The Interaction of Radiation and Matter: 
Semiclassical Theory 



I. REVIEW OF BASIC QUANTUM MECHANICS: CONCEPTS, POSTULATES AND 
NOTATION : 

At the outset, let us, briefly, reconsider why quantum mechanics is necessary? 

• Forces known in classical electrodj^amics cannot account for the remarkable 
stability of atoms and molecules. 

• Disturbed dynamic systems radiate only frequencies which may be expressed as 
differences between certain values (Ritz's Combination Law of Spectroscopy).! 

• All physical systems — viz- material "particles" and electromagnetic fields — exhibit 
wave-particle duality — i.e. to explain particular observations the system must in 
some instances be modeled as a particle and in others as a wave. 

• There is a limit below which the disturbance associated with a observation is not 
negligible and, thus, there is an unavoidable indeterminacy in the prediction of 
observed results. 

To proceed, we recall an ancient comment of P. A. M. Dirac:^ 

"Quantum mechanics. . .requires the states of a dynamic system and the dynamical 
variables to be to be interconnected in quite strange ways that are unintelligible from 
the classical standpoint. The states and dynamic variables have to be represented by 
mathematical quantities of different natures from those ordinarily used in physics" 

Following the Master we begin with the general quantum mechanical principle of 
superposition of states. To quote him once more: 



According to classical theory, a disturbed system should radiate certain fundamental frequencies and their 

harmonics. Each fundamental frequency would be associated with one of the systems degrees of freedom. 

P. A. M. Dirac, The Principles of Quantum Mechanics (Revised fourth edition), Oxford University Press (1967). 
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"The non-classical nature of the superposition process is brought out clearly if we 
consider the superposition of two states, A and B, such that there exists an 
observation which, when made on the system in state A, is certain to lead to one 
particular result, a say, and when made on the system in state B is certain to lead to 
some different result, b say. What will be the result of the observation when made 
on the system in the superposed state? The answer is that the result will be 
sometimes a and sometimes b, according to a probability law depending on the 
relative weights of A and B in the superposition process. It will never be different 
from both a and b. The intermediate character of the state formed by 
superposition thus expresses itself through the probability of a particular result 
for an observation being intermediate between the corresponding probabilities 
for the original states, not through the result itself being intermediate between 
the corresponding results for the original states. 

"...The superposition process is a kind of additive process and implies that states 
can in some way be added to give new states. The states must therefore be 
connected with mathematical quantities of a kind which can be added together to 
give other quantities of the same kind. The most obvious of such quantities are 
vectors" 

With this motivation, Dirac introduced the whimsical name ket and the ingenious right 
half-bracket notation - I ) - to represent the vectors connected to the states of a quantum 
mechanical system.^ Thus, the superposition of two states to form a third is represented in 
abstract ket vector notation, as 

\R)=cM) + c^B) [I-la] 
where the c„ 's may be complex numbers. More generally. 



The state is specified by the direction of a ket vector and any length assigned to the vector is irrelevant. 

R. Victor Jones, March 6, 2000 
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For any set of vectors one can always define a complementary set of dual vectors. In 
Dirac's development, the dual vectors are called bra vectors and denoted by the other half 
of the bracket symbol |\ l| . A bra vector {P\ is completely defined when its scalar product 

— i.e. \P\\Q) = \P\Q) — with every ket vector 1 2/ is known. The scalar product is linear 
(or anti-linear) in the sense that: 

{B\{\A)+\A')} = {B\A) + {B\A') [I-2a] 

{B\{c\A)}=c{B\A) [I-2b] 

{{B\ +{B'\} I A) = {B\ A) + {B' I A) [ I-2c ] 

{c(B\}\A) = c{B\A) [I-2d] 

If (P\A) = for all |A), then (P\ = 0. [ I-2e ] 

Although bras and kets are very different creatures, it is assumed that there is a one-to- 
one correspondence between bras and kets such that the bra corresponding to the ket 
|A)-i-|A') is the sum of the bras corresponding to the kets |A) and |A') -- viz. (A|-i-(A'| 

— and the bra corresponding to the ket c|A/ is the product of c* times the bra 
corresponding to the ket |A) — viz. c*(A|.4 It is also assumed (an especially 



Dirac uses the words conjugate imaginary to connote a corresponding bra-ket pair and, thereby, draws 
attention to an important distinction. The words conjugate complex are reserved for the pairing of ordinary 
complex quantities which can be split into real and imaginary parts. In particular, the real part of a complex 
quantity is given by one half of the sum of the quantity and its conjugate complex. However, bras and kets are 
representations in dual vector spaces and cannot be added! 



R. Victor Jones, March 6, 2000 
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important assumption) that \P\Q) = \Q\P) - so that \P\P/ is real ~ and that 
(P|P>>0.5 

To briefly recapitulate, the directions of vectors of a dual abstract vector space are to be 
associated with states of a dynamic system. To extend the Diracian view, dynamic 
variables and observables of the system are to be associated with (real) abstract linear 
operators which transform one vector into an other ~ viz. 

\F) = Op\A). [1-3] 

Linearity requires that Op{\A) + \A')} ^Op \A)+Op \A') [I-4a] 

and Op{c \A)} =cOp\A) [ 1.4b ] 

The conjugate imaginary or adjoint of \F) is obtained by operating (to the left) with the 
hadjoint of the linear operator ^ on the conjugate imaginary of | A) ~ viz. 

{F\ = {\F)f = {Op\Af = {A\Op\ [1.5] 

By assumption \P\Q) = \Q\Pf so that 

{A\0p^\I^,= {F\0p\A)* [1-6] 
From which we may easily show that 

Op "'""^= Op [ I-7a ] 



^ We may thus identify the numbers {P\2 ) as elements of a hermitian matrix - Charles Hermite (1822-1901). 

^ If a hnear operator corresponds to a dynamic variable, then the adjoint of that operator corresponds to the 
conjugate complex of the dynamic variable. 

R. Victor Jones, March 6, 2000 
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and (op I Op,) =0p. Op, [ I.yb ] 

To further develop the notion of linear operators, we examine the eigenvalue equation 

Jl\F)=a\F) [1-8] 

where is a Unear operator and a is a number. In most instances, the equation arises in 
a context wherein is a known linear operator, which corresponds to a dynamic 
variable, and the number a — an eigenvalue — and the ket \F) — an eigenvector -- are 
unknowns to be determined. Only self-adjoint (i.e. real or Hermitian) linear operators are 
of any use quantum mechanics so we may write quite generally 

Jl\P) = p\P) [I-9a] 

{Q\J^=q{Q\ [I-9b] 

It follows that all eigenvalues are real numbers, that eigenvalues associated with eigenkets 
are the same as those associated with eigenbras, and that conjugate imaginary of any 
eigenket is an eigenbra belonging to the same eigenvalue (and conversely). If we may 
assume that a complete set of eigenvalues and eigenvectors of a given dynamic variable ^ 
is known, we may use the following exquisitely self-explanatory notation to delineate the 
manifold of all eigenvalue equations: 

m)=m)-^ ^io=no; mi=nn-^ [i-ioi 

It is easy to show that within this manifold, the eigenvectors belonging to 
different eigenvalues are orthogonal! 

Representations and Quantum Conditions ; 

General notions : 

R. Victor Jones, March 6, 2000 
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As pointed out earlier, a bra vector is completely defined in terms of a set of numbers ~ 
the representatives -- which are the values of the scalar product of the bra with 
every ket in a representation (and conversely). Thus, if we have a complete set of 
discrete, orthonormal state vectors |t|) we may write a superposition of discrete states 

as 

11 

If we have a complete set of continuous, orthonormal state ket vectors |r|) we may 
write a superposition of continuous states as 

Consideration of these two expansions leads to extremely valuable expressions for 
expansions of the identity operator ~ viz. 

^=Lll)<Tll [I-12a] 

11 



or 



This is, essentially, a diagonal dyadic representation of the identity operator. More 
generally, a linear operator has the following dyadic representation 

Op=II|^)<a<^10p|0 [I-13a] 

Op = JJ i^-xrKriop [I-13b] 



R. Victor Jones, March 6, 2000 
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Since commuting linear operators have simultaneous eigenstates, it is a critical matter to 
know whether or not two observables commute and, if they do not, to know the 
commutation relation needed to replace the classical commutative law of 
multiplication. Fortunately, the so called method of classical analogy provides a 
systematic means for establishing such relationships for the very large class of dynamic 
systems which obey the Euler-Lagrange equations of motion.^ If such a system has n 
degrees of freedom, it is describable in terms of n pairs of canonically conjugate 
coordinates and momenta p^ . ^ The essence of the method of classical analogy is 
embedded in statement that the commutator is the quantum mechanical 
equivalent of the classical Poisson Bracket and is given by 



Thus, the fundamental quantum conditions or commutation relationships are 



This class includes fields, in particular the electromagnetic field, as well as the mechanics of particle motion. 
Precisely: A pair of dynamic variables r and s are said to be canonically conjugate if there is function of the 




[1-14] 



[I-15a] 



two variables such that 



dr 



ds 

and — = — 
dt 



dt ds 



dr ' 



For a second function of the canonical variables we have 



JiF d^F dr d^F ds 3^ dfH 



dt dr dt ds dt dr ds ds dr 



PB 



where {j^ , !B }pg is the classical Poisson Bracket. 
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[Pr^Ps] = PrP.-P.Pr=0 [ I-15b ] 

[?r. Ps]=^rPs-Ps^r=i^ {lr ' Ps }pB = ^ [ I" 15c ] 

These relationships provide the foundation for the analogy between quantum and 
classical mechanics and they show that classical mechanics may be regarded as the 
limiting case of quantum mechanics when fi tends to zero. 

The Schrodinger or Coordinate representation : 

The canonical coordinate(s) is (are) diagonal in the Schrodinger representation so that 

1/ // /// \ ////// I / // /// \ r T 1 T 

?1 ?2 ?3 •/= ?1?2?3 m ?2?3 / [1-16] 

For simplicity, consider a system with a single degree of freedom. For such system 
any ket has a representative ^q') = ~ i-e. 

W{q))= J dq\q')(q\yv{q))= J dq\q)yv{q) [ ] 

The linear operator d/dq acting on this ket generates a new ket which has a 
representative which is the derivative of the first representative ~ i.e. 



d\\t{q)\ 

— T^) [I-18a] 
dq I 



which means that 



R. Victor Jones, March 6, 2000 
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—j=jdq\q')[q 



T I dq 



dq I 



[I-18b] 



Now the crucial commutator of q and djdq may be obtain by the following argument: 



f^«W«))=|l«*(«))= 



dqM^q) \_ 
dq I 



[I-19a] 



I . d d 

and, since \m q)) is arbitrary, — q- q— = 1 
' ^ " dq dq 



[ I-19b ] 



Therefore, to be consistent with Equation [ I- 15c ] we take 



p^ = -in — 



[1-20] 



THE MOMENTUM REPRESENTATION : 

The canonical momenta are diagonal in the momentum representationi^ so that 



PiPi Ps 



\PyPi Pi 



// /// 



7 = P\ Pi P3 



1/ // /// 
Pi plpi- 



[1-21] 



By using Equations [ I-18b ] and [ 1-20 ], we can show, for one degree of freedom, 
that the representative of the momentum in the Schrodinger representation is given by 



\ d \ _ d d 

which means that i~ f = ~ and that — is apurely 

[dqj dq dq 

imaginary linear operator. 

1^ The momentum representation is, generally, less useful than the Schrodinger representation. Because of the 
manifest symmetry between q and p we can now easily show that q^= ih djdp^ . 



I f d I d^[q) 
^ It is easily shown that — ', — 

' dq \ dq 
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{q\p\p')= p'{q\p')=-in{q\^^p')=-in ^A^'\p') [ 1-22 ] 

which has the solution \q\p') = c{p') exp ( i p' q / ^) . [ 1-23 ] 

This connection between the Schrodinger and momentum representations leads to the 
important insight: the coordinate (or momentum) representative of some arbitrary ket 
I A) is given by the Fourier components of the momentum (or coordinate) representative 
— i.e. 

'\A) = (27t ^)"^ J exp(-i p'q/Ti) dq {q\A) [ I-24a ] 



{P 



{q\A) = (271 /z)"^ J "exp (t p'q/fl^ dp' {p'\A) [ I-24b ] 

Application of this transform pair to a wave packets in the Schrodinger representation 
leads directly to the Heisenberg Uncertainty Principle. 

A NOTE ON UNITARY TRANSFORMATIONS 

Let 11 be any linear operator that has a reciprocal 11 ~\ Consider the transformation 

a^ = u au-' [1.25] 



where cx is an arbitrary linear operator. This may be interpreted as expressing a 

T 

transformation from any linear operator OC to a corresponding operator OC . The following 
properties hold: 

a. cc^ has the same eigenvalues as the corresponding a ~ i.e. if a\(i)= (i\a') then 
ii a|a')=ii au~^u\a)=ay^ 'a\a)= a'a\a) [l-26a] 



R. Victor Jones, March 6, 2000 
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b. The fundamental algebraic processes of addition and multiplication are left invariant by 
the transformation ~ e.g., 

(ttj a2y = u tti a^u ~^=u a^u a^u a^^ a^^ ^ ^'^^^ -' 

c. If U transforms a real a into a real 

then a^u ='U a ^ 11* = all * 

and u* a^tl = ll*u a =^ u*a^u = au*u 

therefore 11*11= I ^ U*=ll~^ [ I-26c ] 

Such a transformation is say to be a unitary transformation . 
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II. REVIEW OF BASIC Quantum Mechanics; Dynamic Behavior of 
Quantum Systems: 

Quantum Mechanical Equations of Motion; 



To this point our review, we have been concerned with describing the states of a 
system at one instant of time. The complete dynamical theory must describe, of course, 
connections between different instants of time. 

"When one makes an observation on the dynamical system, the state of the system 
gets changed in an unpredictable way, but in between observations causality 
applies, in quantum mechanics as in classical mechanics, and the system is 
governed by equations of motion which make the state at one time determined the 
state at a later time." i° 

Thus, it is only the disturbance caused by the interaction of a system with a measuring 
device that makes the system's behavior cease to be strictly causal. 

Schrodinger Equation of Motion: 

Consider the time evolution of a particular state of an undisturbed system. To deal with 
such a dynamical system, we need a linear operator of the formal 



From Section 27 of P. A. M. Dirac, The Principles of Quantum Mechanics (Revised fourth edition), Oxford 

University Press (1967). 

This is first member of a class of "displacement " operators that we may treat in a similar fashion. See Section 
25 of P. A. M. Dirac, The Principles of Quantum Mechanics (Revised fourth edition), Oxford University Press 




[II-l] 



(1967). 
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Passing to the limit t^t^, this operator yields a related linear operator for the time 
derivative of the state vector with respect to 

|;l^(0)= {^i^^Ti^ Op (0 [ n-2a] 

If it is postulated that Op{t) = [ifi) ^ ^ [t) where ^ [t) is the total Hamiltonian 

(energy) of the system,i2 we obtain the Schrodinger equation of motion in abstract form 
-- viz. 

^^^M0)=^(0|¥(0) [n-2b] 

Equations [ II-l ] and [ H-lb ] are consistent if 

ih-^/T{t,t,)=fK{t)^{t,t,) . [II-3] 
at 

In the Schrodinger representation, the abstract Schrodinger equation of 
motion becomes^^ 

ih -\\f(?,t) = {7\!H (?) |;) r) [ II-4 ] 

Heisenberg Equation of Motion: 

In the Schrodinger picture, as outlined above, we picture the states of the 
undisturbed motion by associating each state with a moving ket, the state at any time 

There are two justifications of this postulate: (a) analogy with classical mechanics (see Equation [ 11-6 ]) and (b) 
relativistic invariance vw-a-vw Equation [ 1-20 ]. 

Since \y^[t))= \\\ dr\r)[r\^^{t))= \\\ dr\r)y^iCr ,t) 

R. Victor Jones, March 7, 2000 
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corresponding to the ket at that time. In the Heisenberg picture a unitary 
transformation is applied which brings to rest the kets corresponding to states of 
undisturbed motion. In this picture, the appropriate equation of motion is one 
describing the motion or evolution of linear operators (dynamic variables) ~ viz. 

Op =1 ~'0p T or T0p, = 0pT [II-5a] 

where Op is a fixed Schrddinger dynamic variable and Op , is a time varying 
Heisenberg dynamic variable. Differentiating with respect to time 

at at at 

and then using Equation [ II-3 ], we obtain the Heisenberg form of the equation of 
motion as 

ih^^^ = T ~'op y{ T op, = [op„y{,] [11-6] 

where = ^ ~^ y{ T . Note that this equation ~ i.e. Equation [ II-6 ] ~ is 
consistent with the classical analogy discussed earlier. 

ENERGY Eigenvector Representation - The Heisenberg Representation : 
In the Heisenberg picture the stationary states l^^) correspond to fixed eigenvectors of 

a time independent Hamiltonian . 

Case 1 ~ Time independent Hamiltonian: 

To study the time evolution of a given state when the Hamiltonian is time 
independent, we expand the state vector in terms of these fixed energy eigenvectors 
and write Equation [ II-2b ] as 
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ihjW{t))=lH |v|/(0) = li^ K)(£„|¥(0>=I^J^«)(^XO) [n-7a] 

n n 

and ih^ '\KHt))=LK' {K\K^KHt))=K (^X?)) • [ n-7b ] 

Therefore 



(£„|v|/(^)) = (£„|v|/(0))exp^ =(^X0))exF(./(O/) ^^^^^ 



[11-9] 



= C„(0)exi^-/a)„?) 

We may expand an arbitrary wave function in terms of these eigenvectors ~ viz. 
X|/(F,?) = (?|v|/(?))= YJ\A^n)'\KH^)= L{^«W0))(F|£„)exp(. mj) 

n n 

= L<^n(0) "«(^) exp(-/a)„?) 

where M„(r) = !£■„) and the coefficients C„(0) = (£'Jv|/(0)) are, in this 
instance, of course independent of time. 

Case 2 ~ Time dependent Hamiltonian: 

Consider now a time dependent Hamiltonian in the form 

H{{) = H ^ + Hi{) [11-10] 
In this instance, we may write Equation [ II-2b ] as 

«^ n n [n-lla] 
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which becomes, on operating from the left with {E^, 

at n n [Il-llb] 

If we write (£'„|\i/(?)) = C„(?) exp(- mj) - see Equation [ II-8 ] - then 

[/^C„(0+£„C„(?)]exp(-/a)„?) 

= £„C„(?)exp(-iQ)„r) + X;(£;|^i(?)|^„')QWexp(-iM 

n 

Therefore, when we expand an arbitrary wave function a la Equation [ II-9 ], we 
find 

{r|vi/(?))=l^(f|£„X^Jxi/(?)) = £c„(?)(?|£jexp(-/(o/)=Ec„(?) M„(?) exp(-ja)/) [ll-12a] 

n n n 

and C„(r) = --^ X^C„,(r) {EX>iii) exp( / co„,r) [ II-12b ] 



where co =0) -co - and 

rm- n n 

'^E^Hit)\E,)= JJJ dru:\f) Hit) .„fr). [IMS] 
If we write cj^t) = C„(?) exp(- j (Oj) , c„{t) obeys the equation of motion 

c„(?)= -/(0„ c„{t)--J^c„{t) {E„\H,{t)\E„) [ 11-14 ] 
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Time Dependent Perturbation Theory; 

Elementary ideas - First order iteration: 

In set of equations denoted as Equation [ II- 12b ] assume that 

C,.(0) = 1 and C„ (0) = for n ^ / 



so that in first-order 



Cp) ^cHXt) =- i {E„\!Hit)\E,) exi( i (0, t) 



[11-15] 



[11-16] 



Suppose that 



... , , fexpf/co ?) + exp(-/co 



[11-17] 



for t >0, then 



2h 1 /(a)„,+G),) 



[11-18] 



Case 1: co =0 



HKKm" sin^(a),V2) 



[11-19] 



which is valid as long as C,.(r) ~ 1 . Equation [ 11-18 ] shows that transitions 
are more likely if energy is conserved between initial and final states. 

Case 2: ca^^O 

In the so called rotating-wave approximation we neglect the first term in 
Equation [ 11-18 ] so that 
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|(^„N^,)f sin^[(a)„,-co„)^/2] 



= ^2 —rr- [ 11-20 ] 

which is again valid as long as C,(?) =1 . This equation shows that transitions 
are unlikely unless the resonance condition co^ = (O^. is satisfied. 

Let us suppose that we have continuum of energy levels represented by a density of 
states function p (co) . Therefore the total probability for transition out of the initial 
state has the approximate value 

r . M (ivM^ f . 2 sin^[(a)-a),)?/2] 

4^ [(to-(0,)V2] 



For (to-to^) t« % 



2 



Py « J p((o) ' Jo) [ II-21b ] 



d Pr 

so that initial the transition rate ~ i.e., — - — is linearly dependent on t. For 

dt 

2 

longer times, it is reasonable to assume that the "frequency-width" of p(co) ^ ((o)| 
is large compared to the inverse of the elapsed time — or more precisely — —we 

may further approximate the total transition probability as 

Pr-p(co.)t^^ ^^«-^p(o).)r^^(co.)|% [n-21c] 
which yields the famous Fermi Golden Rule for the transition rate 
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Higher-order Time Dependent Perturbation Theory : 

Consider once again perturbation solutions for a Hamiltonian of the form 
9{ [i) = 0-C Q + !H -l^i) . At outset we transform our equations to the halfway land of 

the interaction picture by first rewriting Equation [ II-2b ] to emphasize that it deals 
with a state vector in the Schrodinger picture ~ viz. 

i^J^\^s{t))=^ \^s{t)) [II-2b'] 

Then we make the unitary transformation 

|¥.(0) = '2'o(^) K(0)=exp(-/i^o tin)\y^,{t)) [11-23] 
so that Equation [ n-2b' ] becomes 

'^^^ "^oit) JfWr{t))= 'T.it) fH:,\y\f,{t)) + fH:,{t) 'r,{t)\^^f,{t)) [ 11-24 ] 

Using a version of Equation [ II-3 ] we obtain the equation of motion for the state in the 
interaction picture ~ viz. 

If we write 

|¥/(?))='3^/(?)|¥/(0)) [II-26a] 
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so that 



We may fonnally integrate this equation to obtain 



[ II-26b ] 



[ II-27a ] 



By iterating once we have 



1,(0 = 1 + 



dt^ !H^{t) I dt^ i^/(?2) "T^Ah) 



or by successive iterations we have to an arbitrary level of precision 



\dt,!H,'{k)]dh!HI{t,)]dt,!H'{h) 



+ ■ 



DENSITY Operator (Matrix); 

In classical theory any state of a dynamic system is represented as a point in the phase 
space whose number of dimensions is twice (coordinate and momentum) the number 
of degrees of freedom in the system. This point will move according to the classical 
equations of motion. If the state of the system is defined by some probabilistic 
specification, then we know only the probability that the system may be assigned a 
given phase point at a particular time. We may envisage the time varying probabilistic 
specification as a/Zu/rf of density p(^,?2'93'""A'P2'P3'""'^) moving through phase 



[ II-27b ] 



[ II-27C ] 
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space. Each particle of the fluid will move according to the system's equations of 
motion. The following conservation relationship holds: 



dt 



+ ■ 



dPj_ 
dt 



■=-{P>^}pB 



[ 11-28 ] 



John von Neumann first introduced a corresponding density function into quantum 
mechanics. Dirac pointed out that the existence of such a quantum mechanical density 
is surprising in view of the fact that phase space has no meaning in quantum mechanics 
since numerical values cannot be assigned simultaneously to the coordinates and 
momenta. The quantum mechanical density operator (dyadic) is defined as 



p = £|m)P^(m| 



[ 11-29 ] 



where is the probability of the system being in the state |m) . 

The equation of motion for the density operator is easily determined from the 
Schrodinger equation ~ i.e. Equation [ n-2b ] ~ as 



m 

n 



m)(m|- |m)(m|5/^ } 



[ 11-30 ] 



This equation is, thus, the classical analogy of Equation [ 11-28 ]. Let us now express 
the expectation value of an operator Op in terms of the density operator ~ viz. 
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{Op) = £p^<\|/|Op |x|/) = L^vL<¥|Cp |mXm|\|/) 

\\l \f m 

= LL^v(H¥)(¥|Op \m)=Y,{m^Op \m) [11-31 ] 

m \^ m 

= tr(pOp) 

Gibbs showed that when a dynamic system is in thermodynamic equilibrium with its 
surroundings at a given temperature T, the density is give by 

p= (constant) exp[- ^ /kT] , [ n-32a ] 

This formula is taken over unchanged into quantum mechanics. 

p = (constant) exp[- iH /kl] [ 11.32b ] 
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Appendix - Time Independent Perturbation theory 

Suppose that we have a Hamiltonian ^ = 'S+S where !B is a big part and 5 is a small part. Of 
course, we want the solution to the complete eigenvalue problem - i.e. 

!H:\fH:") = fH:"\fH:") [A-1] 

To obtain an approximate solution we expand the eigenvalues and eigenvectors in a series of terms 
of increasing smallness - viz. 

=J{l + ^Kl+^; + ^;+-- [A-2a] 
\9<") = \tH; ) + ) + ) + ••• [A-2b] 
{t{''\ = {j{; \ + \ + ! + ••• [A-2c] 
Let us first consider normalization - i.e. 

(i^"|i^") = ((i^o" \+\^"\+\^2 l + ---){K )+K"> + k2 ) + •••) = ! [ A-3a] 
Compare terms of equal smallness so that 

{fH:"\^") = (i//; l^r; ) =i [ A-3b ] 

(i^o l^r ) +\^r l< ) =0 [ A-3c ] 

{< ) + I < ) + {^2 ko ) = etc. [ A-3d ] 

From Equation [ A-3c ] we see that we may take (i^j ) = without introducing error. Next 
we write an expanded version of the eigenvalue equation [ A-1 ] 
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(!B+5)(K)+K)+K> + -) 

= (i^o"+5/; +5/2" + ^3 +-"){Ko )+Kr z+K" )+•••) 

Again compare tenns of equal smallness so that 

^\-^o /- -^0 1-^0 / 



[ A-4a ] 



[ A-4b ] 
[ A-4c ] 
[ A-4d ] 



Operate through Equation [ A-4c ] with {p{l \ and using Equations [ A- 3b ] and [ A-4b ] we find 

Operate through Equation [ A-4c ] with (i^™ | and using again Equations [ A- 3b ] and [ A-4b ] we 
find 



or we have the representative 



[ A-6a ] 



[ A-6b ] 



and 



\nrn\_y\nr^ Mnl'n |^« \_ \-| \ ) 

m mtn / 



[ A-6c ] 



Operate through Equation [ A-4d ] with {p{l \ and, of course, using Equations [ A-3b ] and 
[ A-4b ] we find 
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or 



yi-o Pl-^i /--^i \-^o 1-^1 / + -^2 --^2 
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[A-7a] 



[A-7b] 
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III. REVIEW OF BASIC Quantum Mechanics: Two- Level 

Quantum Systems: 

The literature of quantum optics and laser spectroscopy abounds with discussions of the 
two-level (two-state) system. This emphasis comes about because the interaction of such 
systems with the electromagnetic field may be treated in great detail to obtain valuable 
analytic results and, hopefully, the analysis of two-level systems generates insights that 
may be extended to more realistic situations. Fortunately, there are several important 
instances in which the application of the two-level model provides a very good 
approximation to a more complete theory. In the following, we label the upper level of the 
system by the letter a and the lower by the letter b. From Equation [ I- 12a ] we write, 
specifically, the wave function of the two level system as 

^r,t) = {r\\^f{t))= C^t) uj^r) exp{-i (Oj) + C,{t) u,(f) exp{-i(0,t) [ HI-l ] 

where we know from Equation [ I-12b ] that the time varying coefficients satisfy, in 
general, the following equations: 

Q(^)=-^{c,(?)(£,|i^i(?)|£,) + Q(?)<£,|i^,(OK>exp( [ni-2b] 

If we take the interaction to be the electric dipole interaction with an applied electric field we 
may write 

^i(0=~Ps"E(Rs>^)=-ef-E(R5,?) [in-3a] 



R. Victor Jones, March 9, 2000 



The Interaction of Radiation and Matter: Semiclassical Theory Page 27 



where denotes the position of the center of the two-level system or atom.!"* Thus we 
(£, \!H it) I E\ =<U,.=-{E,\ei\E\-t(^,,t] [ ni-3b ] 



write 



In all but the most bizarre circumstances we may use persuasive symmetry arguments to 
reason that 

(£,|e?|£,>0 

Thus Equations [ 111-2 ] reduce to 

Cit) =- -^C^t) exp( mj) [ III-4a ] 

Clt) =- -jCXt) '^a^ exp(-/a),,?) [ ni-4b ] 

where V =-(£je x |^,)- e(r, =- p e(r, . 

rabi Flopping -- without damping; 

For an oscillatory applied field 

'^ah =-|pEjcosco^r=-^pE„exp(-ico^r) + c.c. [III-5 ] 

we see, in the rotating-wave approximation, that 



<^iO =^^E„C,(?) exp[ / (a),,-0)J?] 
= ^^0 C,(?)exp[/(a),,-a),)?] 



[ in-6a ] 



The use of this form of interaction needs considerable elaboration, but we defer that discussion until later. 
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An 

= \^fcXt) exp[-/(a),,-a),)?] 



[ ni-6b ] 



where Q^ = p E^/^ defines the so called Rabi flopping frequency which is, of 
course, a measure of the strength of the electromagnetic interaction. Clearly, the coupling 
terms have maximum effect when the frequency of the applied field is resonant with the 
level splitting. In most treatments the frequency detuning of the field is expressed as 
6a) = co^j -co^ and the system's wave function - i.e. Equation [ ni-1 ] - is written in terms 

of slightly modified time varying coefficients by transforming to the rotating frame of 
reference - viz. 




[ in-7a ] 




[ ni-7b ] 



,f) = CXt) exp[/ (-^Sco-CO J t\ ulr) + CXt) exp[/ (--^Sco-CO J t\ u^r) [ m-Vc ] 



The coupling terms in this rotating frame of reference become 




[ in-8a ] 




[ lll-8b ] 



or in matrix form 




[ III-8C ] 
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We look for a solution in the form C{t) = C{0) exp{^ I iCl^ tj where Q.^ is the 
generalization of the Rabi flopping frequency. Therefore, the condition 

det[M-^2'*I] = 



yields the generalized Rabi flopping frequency 

and the general non-damped time evolving wave function 



sm 



cos(in^) + ^sin{ir2^)j 



A(o) 



[ III-9a ] 



[ III-9b ] 



Rabi Flopping -- with Damping : 

Neglected interactions (e.g. spontaneous emission, collisions, thermal fluctuations) limit 
lifetime of a state of a two-level system. One class of lifetime limiting interactions may be 
described phenomenologically by adding decay terms to the equations of motion - i.e. 
to Equations [ III- 8 ] - as follows: 



or 



dt 



d 




( 
i ' 


dt 


(CM) 





-(5(0- /yJ 

ilf (5(0+ iY.) 



A(0 



[ m-lOa ] 

[ m-iob ] 



= ^M'C{t) [in-ioc] 
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Again we look for a solution in the fonn C[t) - C(0) exp( \ i At) so that 

det[M-Al] = 



or 



Therefore 



[in-lla] 



~ where the we refer to Yafe= ilYa"*" yJ average decay rate constant and 



n'*= J[6a)-it(Y«-Y.)]'+t^^| 



as the generalized complex Rabi flopping 

frequency — and the general time evolving wave function — in the rotation frame — may 
be written 



Ait) J 



cos 



(iU^t)- 



sm 



(lU^t) 



COS 



(lU-t) 



+ 



^si„(in«,) 

KY„-y„) + '8ra 



1 



sm 



(iU^t) 



exp(--2Yai t) 



'cXo) 

^Q(O) 



[m-iib] 



Density Matrix Treatment of a Two-Level Systems: 



Recall from Equation [ 11-30 ] 



n n 



Using |aX«l + I^X^I = 1 we may write 



P = -^{p(\a){a\ + \b){b\)^ -H (|a)(a|+|l;)(Z7|)p} [in-12a] 
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If 9{ = , then 

p=i{p[£ja)(a|+£,|&X^| + |aX«Ki + |^)(^K.] 

-[E^\a){a\ + E,\b){b\+fH:,\a){a\ + ^,\bXb\]p} 



Writing matrix elements ~ i.e. representatives ~ of the density operator 

{a\p\a)=-j,{{a\p\b){b\f^C,\a)-{a\i^,\b){b\p\a)} [IH-lSa] 

(Z7|p|Z7>=-^{(Z7|p|a)(a|itfJZ.>-<Z.|5/JaXa|p|^)} [ III- 13b ] 



(«|p|^>=-|{-(^«-^.)(«|p|^) + [(«|p|«)-(^|p|^)]<«Ki|^) 

+ {a\p\b)[(b\^ ,\b)-{a\^ ,\a)]} 



We may again assume, by symmetry arguments, that {a\ \a) = \b\ ^^\b)=0 so that the 
equations of motion for the elements of the density reduce to 

P..-P««= M"'2l^,,p,, + c.c.] [in-14a] 

Pab=-i(^abPab + i ab[p aa ' P bb] [ HI-Mb ] 

To introduce an element of reality, we add to these equations a pair of the most 
intuitively satisfying damping terms (a useful attribute of density matrix 
formulations) ~ viz. 

Pbb-ppaa=lifi~'^'^abPba+C-c]-^[(pbb-Paa)-(pbb-Paa)^] [HMSa] 



R. Victor Jones, March 9, 2000 



The Interaction of Radiation and Matter: Semiclassical Theory 



Page 32 



[ ni-15b ] 



Transform these equations to a rotating frame by taking Pab~ Pab ^^P (~ ' ^) where 
is assumed to be a slowly varying function of time which satisfies the equations of 

motion 

p,, = -/6a)p,, + /rV^,exp( i(o^t) [p^-pj-yp,, [in-16a] 



Pbb- Paa = [i ^'^ab ^^r ^ P ba + ^c] 



-r 



{pbb-Paa)-(pbb-Paa)^] 



[ III-16b ] 



£xp (— ico^?) + c.c.j and ignoring terms proportional to exp(+ / 2(0 
-- i.e. the rotating wave approximation - we find that 

P,, = -/(5co-/Y)p,,-/(/2n^) [p„,-p,,] [in-17a] 

P - P «« = h' P + c-c] - r [(p - p - (p - p J [ iii-i7b ] 

Steady state behavior ; If we take all time derivatives in these equations equal to zero, 
we obtain 



Pat 



i^o [p - P _ ( 5co+ / y) [p - P 



[III-18a] 



r 



(Pbb-Paa)=(pbb-Paa). 



5a)'+Y' 



6a)'+Yll+f^„l7rY 



[ ni-18b ] 



or 
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-I 1 1 1 1 1 1 1 1 1 1 

2. 4. 6. S. 10. 

We have graphed Equation [ III- 18b ] at resonance - i.e. 5co= - and, from the following 
expression, we see how the oscillatory polarization is saturated at high electromagnetic 
powers. In general, the macroscopic polarization is then given byi^ 



That is to say 

\P> = (P P)aa +(P P).^, = ^aaPaa+ Kb Pba + Ka P ab + Kb Pbb = Kb Pba + Ka Pab 

since P^^ = P^^ = by symmetry. 
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If we graph Equation [ IE- 18b ] at resonance - i.e. 5a)= - we see how the oscillatory 
polarization is saturated at high electromagnetic powers - viz 

P = N (dipole moment of system) = A/^ p + c.c] [ III- 19a ] 

P= -T—^ — exp(jco, r) + c.c. [Ill- 19b] 

2(6a) +Y j 

N\^'[d(o-iy][p,,-pJ 
P= F T — ^E„exp(j(0, ?) + c.c, [111-19C] 

2^l5a)^+Y^[l+f2„t/rYj 
The Vector Model of the Density Matrix: 



There is a set of arguments by analogy which is exceedingly valuable in treating 
transient excitation problems in optics. The basis for the analogy lies in the fact that 
Equations [ 111-16] are identical in form to the famous Bloch equations of magnetic 
resonance. 1^ If we make a transformation to a rotating frame ~ i.e., 
M+ = M+ exp (+i co^ t) — the Bloch equations take on the form 



1^ Recall from the theory of magnetic resonance 

dM/dt= Y„,,MxH 

where Y^^„ is the gyromagnetic ratio. This equation of motion is greatly simplified if it is written in 
terms of the circular polarizations M+ = + i'M^ and H+ = + i — viz- 

Mi = +^-Y„.,M,H,±/Y™sM,H, 
and M^=i^ [M^H_ - M_Hj 

To these Bloch added the phenomenological longitudinal (thermal) relaxation time 7] and transverse 
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= +/ 5a) ± / Y„>ag exp (±/ (o^t)H^-U, /T^ [ in-20a ] 



M = 



H_ exp(-/ CO, t) + c.c. 



{M,-Mj/r, [ni-20b] 



where 5co =(0 , - y^.^^ . By comparing these Bloch equations with Equations [ III-16 ] 
and [ III- 17 ] we see that we have identical problems if we make the identifications 



i/T, = r 

l/T, = y 

In other words we have the equivalent equation of motion 

M= Mxh^ff- Y(M-zMj-rzM^ 

where the e^/ec^ive Rabi precession field ^2^^^ is given by 



^^ff= x-5ol) z 



[ III-21a ] 
[ ni-21b ] 
[ in-21c ] 
[ III-21d ] 
[ in-21e ] 
[III-21f ] 

[ III-22 ] 



[ ni-23 ] 



dephasing time T2 so that 



M, = +i y^ yuw+i y^^ H, - M,/r, 
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Initial state of system 



X 




'Rabi field" 



Effect of transient excitation 
- precession of state vector 




'Rabi field" 



State after 90° pulse 




"Rabi field" off 
-start of "free decay" period 
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IV. REVIEW OF BASIC Quantum Mechanics; The Thermal 
RADIATION Field 

Planck's radiation law for thermal Sources; 1^ 

To set the stage for subsequent discussions of laser physics and quantization of the 
electromagnetic field we briefly explore the earliest, seminal notions in the quantum theory 
of Ught. As we all remember, in 1900 Planck found that he could account for the measured 
spectral distribution of radiation from a thermal source by postulating that the energies of a 
certain se^ of harmonic oscillators are quantized! Drawing on Planck's success, 
Einstein in 1905 was able to show that the extraordinary features of photoelectric effect 
could be explained by hypothesizing the corpus cularity of the electromagnetic field. 
The crowning triumph of early quantum optics is Einstein's amazingly simple, 
phenomenological theory of 1917 which provided a quantitative basis for analyzing the 
absorption and emission of light by atoms. 

In treating thermal sources the basic assumption is that emitted radiation is a sample of the 
total electromagnetic field ~ viz. all of the modes of a resonator — in thermal equilibrium 
with its material environment — viz. the walls of the resonator. In the traditional treatment 
of the theory of black-body radiation, a particular thermodynamic system is assumed - 
- viz., a hollow resonator which is a cube of length L with has perfectly conducting 
walls. 1'' Thus, to satisfy boundary conditions at the walls, the electric field associated with 
a particular mode of the cavity is given by 



1^ This section and some parts of the following section draw heavily upon discussions in Rodney Loudon's The 
Quantum Theory of Light (2nd edition), Oxford (1983). 

1^ The detail results of black-body radiation theory are indeed sensitive to the assumed character of the 
thermodynamic systemas will be demonstration presently. 
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K" ' (?' t) = E° (k,^ ,J COS {n V, x/L) sin (ji y/l) sin (tc z/l) exp (/ (O,^,^,^ r) 

'^(f, t) = (fc,^, ,J sin (71 V, x/L) cos (ti y/l) sin (71 v, z/l) exp (/ to,^,^,^ r) [ IV-1 ] 
'-(f, t) = E^ [k,^ ,J sin (7t V, x/L) sin (tc y/l) cos (tc z/l) exp (/ co,^, t) 

where {v^, Vy,v^} is a set of positive integers. To satisfy the homogeneous Helmhotz 
equation we must have 

<v,v.= vjL)%(7U vjLy + (n vjLf= (to,^,^,^ /cf [ iV-2 ] 



The cycle-averaged value of the stored energy density associated with the particular 
mode is given by 



dV 



2V 



cavity 



[IV-3] 



For a thermal source, the most significant experimentally measurable object is the noise 
spectrum ~ i.e., the frequency distribution of the stored energy density. To obtain this 
distribution, we take the energy density in the frequency range between co and (0+ rfco — 
viz. 

W[(o) d(0= ^ y „J <=|with ^^^y between CO /c and (a)+ d(o) /c| 



_/ 



X 



Number of modes with frequencies 
between to andto+^/to 



[IV-4] 
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Thus the density of states is defined as 

p^dk = 



Number of modes with frequencies 
between CO and co+ dd) 



[ IV-5a ] 



In the frequency range where X«2L, the following -space argument holds: 

[two polarization~|rFraction of shelll[ Volume in shell - ] 

pjj. dk ~ I _ II li = Number of k states within shell 

[states perk state JLwith vaid states JL Volume of state J 



»2x 



P 



[ IV-5b ] 



This extremely important density-of-states 
construction may be visualized most 
elegantly in two dimension, 




Cut thiDTigh shell in t-space. 



For use in later discussions, we note here that by this same argument the 2D density-of-states is found to be 
p^(2Z)) = Aco/ 7t and the ID value pJ,lD) = L/li c. 
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Following the traditional (Rayleigh- Jeans) argument, we identify a resonator mode as a 
harmonic oscillatori^ and take its average energy to be the classical thermodynamic value 
for a system with two degrees of freedom ~ i.e. (wlk^ v v )/ = V Thus, we 

obtain the famous Rayleigh-Jeans radiation law (which impels us to worry about an 
ultra-violet catastrophe .) 



Wt(co) d(ii = kT 



3 ,^2 
C 7t 



[IV-6] 



Planck's quantization hypothesis; 

Following Planck we set 



1 ^ 

^ (W)= p(0 = '^n/V= [h(ii{n + i)]/v [ IV-7 ] 



where n = 0, 1, 2, Assuming that the resonator mode is in thermal equilibrium 

with its environment, we may use the Boltzmann probability factor 

exp(-|3 'Ej = txp{-'Ejk^T) 



to find the probability P„ that the mode (read oscillator) is thermally excited to an 
energy <£„ 



P. = 



exp(-|3 £„) exp(-|3 n/jO)) 



" l^exp(-P'E„) 2^exp(-|3n^co) 



= exp(-pn^a)){l-exp(-p ^co)} [IV-8] 



The thermal mean value of n is obtained in a similar manner 20 — viz. 



19 



20 



We justify this assertion in detail later: it suffices to note here, in anticipation of that later discussion, that a 
harmonic oscillator and an electromagnetic mode have analogous Hamiltonians. 

/ A \ // \2 / \ 

It is also easy to show that thermal mean-square fluctuation in n is given by \An/j = \n/j +\n/j . 
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exp(-P n^co) 

n 

= ^ co)} In ^Ilexp(-|3 n ^co) j 

= ^ (o)} In ^{l- exp(-p ^co)}"'j 

_ 1 
exp(P ^co)-! 

which may be plotted 
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[IV-9] 



i 
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Equation [ rV-9 ] taken in conjunction with Equations [ rV-5b ] and [ IV-7 ] leads then 
directly to the Planck radiation law 

WT.{a))da) = {n + {)^ 



lermal 



h 0) Jo) 

3 „2 

c n 



1 

+ - 



exp(|3^co)-l 2 



d(ii 



[IV-10] 



From this expression we see that the Rayleigh- Jeans radiation law is correct in the limit 
k^T » fi(0 , but, to our great relief, Planck has staved off the ultraviolet catastrophe. 
We may then plot this famous and an exceedingly important result. (From the 
particular point view of laser physics, the Planck radiation law gives a measure of the 
background radiative thermal noise spectrum.) 
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Einstein's Phenomenological Theory of Radiative Processes -- 
Einstein's "A" and "B" coefficients: 

We return to the two-level system and, with Einstein, make some physically reasonable 
(early quantum mechanical) postulates concern the absorption and emission of light. 
Consider the two-level system shown below. Einstein's postulates^i are embodied in a 
set of rate equations for the level populations- viz. 



Using our semiclassical theory of matter-radiation interaction, we can easily show that = j > 

but, of course, Einstein did not have that well-developed theory at hand. While the stimulated emission and 
absorption processes are reasonably intuitive, it is the introduction of the spontaneous emission process which 
flows from a profound understanding of the processes by which thermal equilibrium is achieved. 
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dNjdt= N^A^ - N, B,^ W{(0^) + B^, W{(0,) 
dNjdt=-N^A^ + N,B,^ W(a)J- N^B„, W(a),) 



[IV-U] 




fj CO 



Spontaneous 
emission 



Stimulated 
absorption 



Stimulated 
emission 



It particularly useful to examine the results of these equations under conditions ~ viz. 
General equilibrium or steady-state condition; 

- B,^ W(co J + B^, W(coJ = [ IV-12a ] 



Thermal equilibrium condition: 



[ IV-12b ] 



Of course, in the latter case we use fact that the level populations A'^ and A^^ at thermal 
equilibrium are related by Boltzmann's law ~ i.e. 

NJN, = [g, exp(-p E,)]/[g^ exp(-|3 £j] = [gjg^] exp(p hoij [ IV-13 ] 



so that 



W,{0i) = A,, [(&/&) exp(|3 hto) B,^ - 5„J 



[IV- 14] 
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Thus, by comparison with Equation [ IV-IO ], we see that 

{gb/ga) Bha=B„b 



and 



^2 3 
K C 



Bab=Kb 



and, by comparison with Equation [ IV-9 ], we see that 



or 



The rate of spontaneous emission A, 1 ^ 

= = T~r = ^'^rP ^ ^ 

The rate of stimulated emission 5^^W^(co) (n)^ 
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[IV- 15a] 
[IV-15b] 

[ IV-16a ] 
[ IV-16b ] 
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V. ELEMENTARY LASER MODELS -- RATE EQUATION APPROACH: 

POPULATION Inversion and "Lasering" in a Three-Level System 

Obviously, the trick in building masers and lasers has been to find means to "invert" the 
equilibrium population of states.22 First conceived by Charles Townes in 1951, Gordon, 
Zeiger and Townes in 1954 demonstrated the first ammonia beam maser. Degenerate 
quantum states of the ammonia molecule (NH3) are split by the tunneling of the nitrogen 
atom across the plane formed by the three hydrogen atoms. In the Columbia maser, 
thermal populations of molecules in the split ground state (23,800 MHz) are spatially 
separated by an electrostatic field to achieve population inversion in the maser resonant 
cavity. While this two-level system was clearly seminal, it was the multi-level paradigm of 



Mirror Flash tube Partially reflecting mirror 




Trigger electrode Quartz tube Ruby rod 



The Configuration of Theodore Maiman's Ruby Laser (1960) 



Recently the subject of lasing has become of considerable interest. See, for Thomas Mossberg's very 
interesting Web Seminar on the subject (http://decryptor.uoregon.edu/~mosswww/moss-seminar.html). 
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Bloembergen's three-level ruby maser which galvanized the search for possible laser 
systems. 

To quote from A Laser Tutorial (http://members.aol.conVWSRNet/tut/tl.htm) 

"Finding substances in which a population inversion can be set up is central to the development of 
new kinds of lasers. The first material used was synthetic ruby. Ruby is crystalline alumina 
(AI2O3) in which a small fraction of the AP* ions have been replaced by chromium ions, Cr'*. It is 
the chromium ions that give rise to the characteristic pink or red colour of ruby and it is in these 
ions that a population inversion is set up in a ruby laser. In a ruby laser, a rod of ruby is irradiated 
with the intense flash of light from xenon-filled flashtubes. Light in the green and blue regions of 
the spectrum is absorbed by chromium ions, raising the energy of electrons of the ions from the 
ground state level to one of the broad F bands of levels. Electrons in the F bands rapidly undergo 
non-radiative transitions to the two metastable E levels. A non-radiative transition does not result 
in the emission of light; the energy released in the transition is dissipated as heat in the ruby 
crystal. The metastable levels are unusual in that they have a relatively long lifetime of about 4 
milliseconds (4 x 10"^ s), the major decay process being a transition from the lower level to the 
ground state. This long lifetime allows a high proportion (more than a half) of the chronaium ions 
to build up in the metastable levels so that a population inversion is set up between these levels 
and the ground state level. This population inversion is the condition required for stimulated 
emission to overcome absorption and so give rise to the amplification of light. In an assembly of 
chromium ions in which a population inversion has been set up, some will decay spontaneously 
to the ground state level enutting red light of wavelength 694.3 imi in the process. This light can 
then interact with other chromium ions that are in the metastable levels causing them to emit light 
of the same wavelength by stimulated emission. As each stimulating photon leads to the emission 
of two photons, the intensity of the light enoitted will build up quickly. This cascade process in 
which photons emitted from excited chromium ions cause stimulated emission from other excited 
ions..." 
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"The raby laser is often referred to as an example of a three-level system. More than three energy 
levels are actually involved but they can be put into three categories. These are; the lower level 
form which pumping takes place, the F levels into which the chromium ions are pumped, and the 
metastable levels from which stimulated emission occurs. Other types of laser operate on a four 
level system and , in general, the mechanism of ampUfication differs for different lasing materials. 
However, in all cases, it is necessary to set up a population inversion so that stimulated emission 
occurs more often than absorption." 

The following rate equation analysis demonstrates the critical conditions necessary to obtain 
sufficient inversion: 

A "Pumped" three-level quantum system 





A 

^31 


A32 ^ 






A 






B2iU((0l) 



E3, N3 



E2, N2 



E,, 



Master rate equations which include "laser" and "pumping" radiation, but in which thermal 
radiation is neglected: 

^N, = N, A,, + N, A,, + {N, -N,) B [/(to,) + (A^, -N,)B,, C/(to J [ V-la ] 
at 



-N,=N, A3, - N, A,, -{N, -N,) B [/(coj 



[V-lb] 
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— N^ ^-N^A^^-N^A^2-{N^-N^)B^^ U{(£)p) [ V-lc] 

dt 

Of course, only two of these equations are independent, since we assume that there are a 
fixed number of systems N distributed over these energy levels. In what follows, we 
define 

R=U{(>ip)B^^ {N,-N,)/N 
as the pumping rate with N^ = N- N^- N2 . 



We now study the steady state solutions of two independent rate equations ~ i.e. 

0=N, A3, - N, A„ - [N, -N,)B,, C/(a) J 

= -N,(A,, + A,,)-{N,-N,)B,, U{(^,)^-N,(A,, + A,,) + RN 

which may be written 

[A32 - 5 21 C/(co J] + [A32 + A21 + 5 21 [/(co J] = N A32 
[^31 + ^32] + [A31 + A32] = (A3, + A32) -RN 



[V-2a] 



[V-2b] 



By messy, but trivial algebra, we find 

^ n[r (A3, + A,i + [/(coj)- (A31 + A3,)(a,i + B,, [/((oJ)] 

[A,,+Aj[A,,+2B,,U{03,)] 

[V-3] 

n[r (A3, ^(co J) + (A3, + A3,)^ ,1 ^(0 J] 



^3, + A32 
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Therefore, the all important population difference may be expressed as 



_ n[r(2A,, + A,,)-A,, (a,,+A,,)] 



or, with regrouping of terms, as 



A21) (Ajj + Ajjjj 



A21 [A31 + A32] 



1 + 



25 21 [/(CO J 



'21 



n[r (2A,,+A 

21) Aji (Ajj + Ajjjj 



A21 [a 31 + A32] 



1 + 



[ V-4a ] 



[ V-4 b] 



where = {c/n) (Aji/ 25 . Since the spatial change in the Poynting vector associated 

with a beam propagating through the laser medium is equal to the rate of energy density lost 
(or gained) by the beam, you can show that 

1 d 



l{z, to) dz 



l{z, co) = -(A^i - N,) [c/n] {h(0B ^jF{oy)/v) 



[V-5] 



where F(co) is the line shape function of the E^- absorption line. Using Equation 
[ V-4b ], we may write 

1 ^(2A32+A2i)-A2i(A3i+A32)/^cAAr^a)5 2iF(co) ^^^^^ 

= r — — I - I TT — [ V-6 ] 



l[z, co) dz 



A2i[A3i+A32j|l + -y— ^ I 



Therefore, we may write 
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dz 



[V-7a] 



where we identify the ^ain of the medium as 

r(2A,, + A,,)-A^,(a,, +A,^)fc\Nh(a B^, F(co) 



^21 [^31 +^32] 



V 



[V-7b] 




o 



Excited atom 



Ground atom 



4lLfo o o o f^a# • o • ^ 



Rear mirror 
(Totally reflecting) 



Stimulated photons 




Ui.SER 



OUTPUT 



Output mirror 
(Partially reflecting) 



Random photons 



A LASER SYSTEM 



Thus, we have a interesting and useful model for the "gain behavior" (growth in intensity 
with position) of a beam as a function of initial intensity - viz. 



In 



^ /(z,co) V /(z,co)-/„ 

.~r~ — 1 — 

V o J crir 



[V-8] 
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For « 4„., l{z,(o) = exp(g z) "Inverse Beer's Law" 

For /„ » l{z, (o) = 4 + /„„ g z Linear regime 

Now consider elementary laser operation with feedback provided by two mirrors separated 
by a distance L. We first calculate the famous formula for threshold gain from the 
following consistency condition: 
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/„ exp(-a 2L) exp{g„ 2L) r,r^=I„ 



[ V-9a ] 



or 



VI '27 



[ V-9b ] 



where a is included to account for residual cavity losses and the 's are the respective 
reflectivity's of the two mirrors. At high power — i.e. at powers well above threshold ~ 
the consistency condition becomes (neglecting residual absorption) 



1 



or 



"1-1 



-1 



[ V-lOa ] 
[ V-lOb ] 



From Equation [ V-7b ] we may write 



where 



A^,(a,, + A,^) ^ „ „ ^ ^21 [^31 + ^32] V 
= — ^ and R,^-R^ = g,y, - 



(2A32 + A21) 
Therefore Equation [ V-lOb ] becomes 



(2A32+A21) Uy A^^co ^21 F(co) 



{ 



1 



= \ r , ^'-f'"^ [{R-Rj = C{R-Rj 

Il'-2 J 



[V-12] 



[ V-Ua] 



[ V-llb] 
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This result, which is plotted below, is archetypal of the threshold behavior all lasers. 

LASER LIGHT OUTPUT "WELL ABOVE" THRESHOLD 



S.T 



6.-- 



4." 



2." 




Optical "pumpi ng" rate 
H 1 1 1 1 1 R 
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Appendix: Common Lasers 
Pink Ruby Laser Laser: the first pulsed laser 

Energy Levels of Cr3+ in a Trigonal Crystal Field 



Broad 
F levels 

Pumping light 





Ground slate 



Rapid 
non-radiative 
M transitions 



: : Metastable levels 



Laser 
transition 




( 694.3 nm ) 



Enerigy levels of chromium ions in ruby 



Absorption Spectrum of Rl and R2 of Cr3+ in Rub 
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HeNe Neutral Gas Laser: the first continuous laser 



160,000 



150,000 



- Is2s 



ls2s 



E 
(J 

LU 



140,000 



130,000 



V 



V 1, 



Is^ 



Odd 
parity 




2f?4p — 



r 



20 



19 



18 ^ 



LU 



17 



16 



2^ ^ 



Mechanism: A dc or rf electrical discharge is established in a gas mixture of helium (= 90%) and 
neon (-10%). The helium atoms are excited to long-lived metastable states by a 
variety of processes involving energetic electron impact. The excitation energy of the 
helium metastables is transfered to neon atoms by resonance scattering to yield laser- 
active neon metastables. 



Particularly strong emissions 



Is^ls^lp^Ss: 'Pj ^ ls^2s^2p^3p:^Si 
ls^2s22p^5s: ^Pj ls22s^2p^3p: '^V^ 



1 



A-vacuo = 0.5435 |xm; 
^vacuo = 0.6330 |xm; 



>.aii= 0.5433 |xm 
>.aii= 0.6328 ^im 
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ls^2s^2p^4s: ^Pj ls^2s^2p^3p: ^vacuo= 1.1526 ^im; X^ = 1.1523 |im 

lsW2p^4s: ^Pj ^ lsW2p^3p: ^Sq A.vacuo= 1-5235 ^im; >.air= 1-5231 |xm 
lsW2p^5s: ^Pj ^ lsW2p^4p:^P2 A-vacuo = 3.3922 ^im; >.aii= 3.3913 |im 

Rare Earth Ion Solid State Lasers: 

Neutral atom electronic configurations: 

noble gas core xeon: [Xe] = ls^2s^2p^3s^3p^3d^^4s^4p^4d^^5s^5p^ 

rare earth atoms: ls^2s^2p^3s^3p^3d^^4sHp^4d^^4f''5s^5p%s^=[Xe] 4f%s^ 

Trivalent ion electronic configurations of common laserable rare earth ions: 



praseodymium (Pr"^^): 


= [Xe]4/2 


Ground state: 


neodymium (Nrf"^): 


= [Xe]4/3 


Ground state: ^19/2 


europium (Eu"^^): 


= [Xe]4/6 


Ground state: 


hohnium (Ho^^): 


= [Xe]4/io 


Ground state: 


erbium (Er^^): 


= [Xe]4/ii 


Ground state: ^li5/2 


thulium (Tm+^): 


= [Xe]4/i2 


Ground state: ^Hg 


ytterbium (Yb"^^): 


= [Xe]4/i3 





Important characteristic of rare earth ions: Outer 5s and 5p shells 
effectively shield inner 4/ electrons from significant interaction with the 
local crystalline field associated with the charges on neighboring ions in 
solid state configurations. 

1. YIG Laser : Nd*^ randomly distributed on Y*^ sites in Y3AI5O12 (garnet) 
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The absortion and emmission cross sections for 
the '*Ii5/2 ^^13/2 transition in Er+^. 
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VI. Semiclassical laser Theory 23 

Laser Self-Consistency Equations 

Lamb's theory of laser operation provides a very powerful means for interpreting and 
predicting complex time dependent behavior without invoking all the intricacy of a full 
quantum mechanical theory. It is semiclassical, self-consistent theory in the following 
sense: 



1 



quantum 
theory 



statistic 
summation 



macroscopic 
Maxwell theory 



self-consistency 



Suppose that the field in the laser is represented in the form 

E{z,t) = \T. £„(0exp[-/{(O„?+(^„(0)]f/„(z) + c.c. [VI-1] 



where, for example, 

in a simple a cavity laser U^[z) °^ sin k^z = sin 



n 71 



and in a ring laser t/„ (z) o« exp ( i z) = exp 



n 271 ^ 



standing wave [ VI-2a ] 



traveling wave [ VI- 2b ] 



An adaptation or interpretation of Lamb's "semiclassical" or "self-consistent" laser theory as first presented in 
Phys. Rev. 134, A1429 (1964) and refined in countless other treatments. In these lecture notes we drawn 
extensively on M. Sargent III, M. O. Scully and W. E. Lamb, Jr., Laser Physics, Addison-Wesley (1974) and 
P. Meystre and M. Sargent III, Elements of Quantum Optics, Springer- Verlag (1992). 
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With this representation for the field the induced polarization can be expressed as 

P{z,t) = ^L^Pn [t] exp[-i(a)„f+(l)„(r))] U,{z) + c.c. [ VI-3 ] 

n 

where fP„ (?) is the complex, slowly varying component of the polarization of the nth 
mode. If we take the wave equation in the form^^ 

ot c ot at 

where second term is included as a means to account for cavity losses. From Equation 
[ VI- 1 ] assuming that £ „ , fP„ , and are slowly varying functions of time^^ 

-V^E ^- = - „ exp (co„ t^K)] 1^ + C.C. 

a Z I „ a z 

^ " [VI-5a] 

= —Y,Ql exp[-/(a)„M;?+(^„)] U^ + c.c. 

n 

where Q.^ = c are the eigenfrequencies of the cold resonator eigenmodes. 

|^^^^7^ = ^L{£„-''£>„+(^„)}exp[-j(03„?+^j]t/„(z) + c.c^ [VI-5b] 



In this formulation we are assuming that V P » . 

The so called slowly-varying amplitude and phase approximation (SVAP) is used extensively in 
treating problems in laser dynamics. In the SVAP approximation it is assumed that 
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+ C.C. 



[ VI-5c ] 



3^ 3^2 = 2^ |'E„-/2 2;„(to„+^„) + '£„[-(to„+(^„) -/^„J|exp[-/(a)„?+^j]c/„(z) + c.c. II-5d] 



which reduces in the S VAP approximation to 

32 E d^E{z,t) 1 



7^ -:;L{-i2£„a)„ + £„[-a)^-2a)„(^j}exp[-/(a)„M;?+(t)„)] C/„4z) + c.c. [VI-5d'] 



and from Equation [ VI-3 ] 

a'p d'p{z,t)_i 



^ = 2^i^''>" {(0„+(^J + ^P„[-{(0„+(^„) - i(i)„J| exp[-i (a)„?+^J] f/„(z) + c.c. 
{-«n ^P„}exp[-/ (o)„?+(^„)] f/„(z) + c.c. 



[VI-6] 



If these representations are to be valid when substituted into Equation [ VI-4 ], the 
following equations must hold: 



n n 



m £ - i 2m £ -co' £ -2m 6„ co' ^„ 

n n n n n n n in n ^ n n 



[VI-7] 



We adjust \hQ fictional conductivity to account for the loss of energy or time decay of 
the given mode ~ viz. we take 0= 8„ co„/g „ where Q n is the "Q" of the mode. Equating 

real and imaginary parts of Equation [ VI-7 ] we obtain Lamb's master laser self- 
consistency equations -- viz. 
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03 


9ot( <? 1 

-'-n I 


r VT 8n 1 
L V l-oa J 


n n n 


£ -2(0 (b £ = 

n n Tn n 




[ VI-8b ] 


or a)„+(^„»r2„ 


2e<, ^ 


}/£„ 


[ VI-8b' ] 



Polarization of the medium ; 

The simplest application of Lamb's theory assumes that the laser medium consists of an 
ensemble of two-level atoms with a single well-defined transition ^(0^;, which is 

homogeneously broaden . Possible "pump" and "decay" processes are indicated 
schematically in the following figure: 




The model describes, in effect, a so called four-level homogeneously broaden laser 
where both levels relaxed to some presumed ground state and are "pumped" via higher 
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energy excited states, p [s,z,t',t) is defined as the density matrix operator at a place z and 

a time t which is associated with atoms excited to the state s ( where s = a or & ) at a time 
f and 'k^[z,t') as the rate at which atoms are excited to the state s (per atom, per second 

per unit volume). The macroscopic polarization is, thus, given by 



P{z,t) = Y,l df l,{z,t)(e mp ) 

sw 

= pY^j dt'\{z,t'}p,^{s,z,t\t) 



[VI-9] 



+ c.c. 



If we project this expression on to the modes of the laser cavity and compare with Equation 
[ VI-3 ], we see that 

fP„(0 = 2 ^exp [ / (a)„ t+,^it))] I dzu: (z) l^jdt'Kizj) p., {s,z,t\ t) [ Vl-lOa ] 



where ^= \ dz \ U^[z) \ is the mode normalization factor. To facilitated the integration 



of this equation, we re-write it in terms of the population matrix operator ^ {z,t) — 



VIZ. 



^it) = 2 pexp [ / (q)„ t+K{t))] ^]dz K (z) (z, t) [ Vl-lOb ] 



n 



Where p {z,t) ^Y^l dt' K{z,Op {s,zy,t) ■ [VI-U] 

S -oo 

The equation of motion of this population matrix is found by differentiating its defining 
equation ~ i.e. 
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=LK{z,t)p{s,zxt)+I,l dt'\{z,Op{s,z,t',t) . [Vl-12] 

s s ^ 

By definition p^p (5,2, =6 „j 6p 3, so that first term on the right-hand-side is replaced by 
the matrix 

X,{z,t)\ 

If we interpret the p {s,z,t',t) components of the second term on the right-hand-side of 
Equation [ VI-12 ], as the time derivatives of the pure state density matrix components 
given by Equations [ III- 14a ] and [ III- 14b ] of this set of lectures notes,^^ the component 
equations of motion for the population matrix p {z,t) become 

Paa = - Ya Paa " '^'ab f^ab + Cc] [ VI-13a ] 

Pbb = ^b - Pbb + [i '^ab Pab + c.c] [ VI-13b ] 

Pab = -('«ab+Yab)pab + '^"'l^ab [paa-Pbb] [ VI-13c ] 

Formal integration of Equation [ VI- 13c ] yields 



Including damping, Equations [ III- 16a ] and [ III- 16b ] of this set of lecture notes would have the form 

Paa = -YaPaa-[^^"''J^abPab + C-C.] 
Pbb = -YbPbb + [i^"''^'abPab + C.C.] 

Pab = 0)ab +Yab) Pab + « ab [p aa " P bb ] 
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Single Mode Operation; 

For a single-mode standing wave field in the rotating- wave approximation 

'^ab(2. =-\p'^Xt) exp [-i {«„ t'+^XtD] U. {z) [ VI-15 ] 

^ / ph-' cxp[-io)j] (z) J '^J/exp|-[/(a),,-(o„)+Y,b](?- t)- i^,{t)^'E ,{t) [paal^-O- Pbb(z>0] 
^ipff' Qxp[-i(xij] U^{z)\^ t/Texpj-[/(a),,-a)„)+YjT-/(^„(?-T)|£„(?-T) [M^'^- '^)] 

The integration can be done simply, if the changes in £ (|) and 1^^,^^ - Pj.^ are negligible 
in a time so that 

Pab(2'^) = -^ '"^^^'^ '^.{t) exp[-/(c0„ „(?))] C/„(Z) [Paa(Z'^)-M2'^)][^"Kb-«„)+Yab] ^ 

^ [VI-17] 
= -- ipr exp[-i(a)„ [paa(^>?)-M^'^)]^^Kb-«n;Yab) 

where 2^(m- viw) = [i (m- v) + w] represents the so called complex Lorentzian 
denominator. Substituting this result into Equations [ VI- 13a ] and [ VI- 13b ] we find 
the previously discussed rate equation approximation for the population components - 
- viz. 
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P'aa= ^a-YaPaa-^(Paa-Pbb) 
f^bb= ^b-YbPbb+^(Paa-Pbb) 



[ VI-18 ] 



with the rate constant is given by 



2 



^^£„^(0|f/„(z)fY.,[K-cojVYL|' 
^ k„^(r)|f/„(z)fY-;,x(a),,-a)„;Y,,) 



[ VI-19 ] 



where X(m -v; w) = [(m -v)^ + w^] represents the dimensionless Lorentzian 
function. 



For steady state ~ i.e. fj^^ = Pbb ~ ^ 



Paa-Pbb 



?^aY:-^bY;' N{z) 



Ya Yb 
lYa+Yby* 



\ + R{z)lRs 



[ VI-20 ] 



where N{z) is the unsaturated population difference and ^5 = Ya Yb (Ya+Yb) 
Combining Equations [ VI- 10b ], [ VI- 17 ] and [ VI-20 ] we obtain 



N{z)\USzt 



R{z)/Rs 



[ VI-21 ] 
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Although this result can be integrated exactly, there is more insight to be gained by 
expanding the safara^/on denominator in powers of R{z)lRs — viz. 



^ii) =- ^"'£„ 2^co,b -«n; Yab) ^]d^ M^) I UXz) f £ 



[ VI-22a ] 



To second order in R[z) /R^ we obtain 



4 " 



Ya+Yb 



V lab y 



4^ab-(0n;Yab) 



[ VI-22b ] 



where 



and 



M')-7jj-J*Mz.')|(^.W 

n 



^^l^^E/ 1 



" 2 Y,Y, 



[ VI-23a ] 



[ VI-23b ] 



is the dimensionless mode intensity. 

To complete the self-consistency loop, we combine our expression for the complex 
mode polarization. Equation [ VI-22b ], and the self-consistency conditions. Equations 
[ VI-8]. 

From Equation [ VI-8a ] 



27 Exactly 



where 



f[w) = - l-(l + w)" 



and w= T 



f \ 

Ya+Yb 



4^ab-(0n;Yab) 
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2 " ne„ 



Ya+Yt 



4^ab-«„;Yabw) 



J 



Multiplying this equation by 9 ^ ^ , we obtain 

lb lb 



I 



CO 

n 



1 ^ 

2 "he^ 



' O - 1" 3 

— Jx(co,,-co„;yJ A^^l--/„ 



^VKb-co„.;Yab)] 
^ Yab ; JJ 



[ VI-24 ] 



[ Vl-25a ] 



This nonlinear equation of motion may be written to advantage as 



[ VI-25b ] 



where the linear gain is measured by 



1 . 



2 "he^ 



V'ab J 



4^ab-(0n;Yab) 



[ Vl-26a ] 



and nonlinear saturation by 



8 ^e,, Jab 



Ya+Yt 



^l0)ab-«„;Yab) 



[ VI-26b ] 



From Equation [ VI-8b' ], we see that 
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Thus 



a)„+(|)„=n„ + -a)„^ ^ , — ^i:(co,b-co„;Yab 



" 2 

= ^„+Cy„-Pn^„ 



Yab 



Ya+Yb 



^(«ab-«„;Yab)| 



[ VI-27 ] 



where the linear mode pulling effect is measured by 



"""^ 2 ^ ITT ~1 — '^K-''-^-) 

^ "*^o lab lab 



[ VI-28a ] 



andnonlinearpMS^/ng effect by 

p _3 a)„p^ Kb-con) 

" 4 ^8„Y,b Yab 



Ya+Yb 



[4^ab-«„;Yab)] 



[ VI-28b ] 



From Equation [ VI-25a ], we find the "at-resonance"threshold condition to be: 



2QJ 2 "M^ 



ab J 



1 p^N, 



Th 



On Ky 



ab 



[ VI-29 ] 



where A^^^, is the required population inversion at threshold. 



For the general steady state condition , we have 



/ = - 



2 r Ya Yb 



- «„/P„ 



[ VI-30a ] 
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or 



4 

I.= - 

" 3 



Yab ^4^ab-»n;Yab)- A^Th/A^ 



Single- MODE, Standing- wave Laser Intensity 



[ VI-30b ] 




The time dependent intensity build-up follows from integrating Equation [ VI-25b ] 
to obtain 

«n /n(0)exp(a„r) 



a„+p„/„(0)[exp(a„f)-l] 



[ VI-29 ] 
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Transient Intensity Buildup 




2 4 6 S 10 

time 



MuLTiMODE Operation ; 

When we consider the possibility that two or more laser modes may be simultaneously 

excited or oscillating, things get quite a bit more complicated and we must refine our 

analysis. In particular, we must treat the time dependence of the population difference 
Paa ~ P^hb with greater care since variations in the population difference tend to modulate 

and couple the possible modal excitations. 

Reflection on the form of Equations [ VI- 13a ] and [ VI- 13b ] in the absence of modal 
excitation, suggests that we take 

N{z,t)^[pT {z,t)-i^ll\z,t)]= X, y-J-K r: [VI-31 ] 
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as the zeroth-order approximation to the time dependence of the population difference. 
With finite excitation of one or more modes we may write in the rotating-wave 
approximation 

'^M'^ =- ^pL2;„(0 exp[-/(a)„ U,{z) [ VI-32 ] 

^ n 

and use Equation [ VI- 14 ] to obtain the following first-order approximation to the time 
dependence of the off-diagonal component of the population matrix: 

= / n-' J j/exp[-(/a),,+Y,,)(^- o]'^a.(^,o[p:H^>^)-p:n^,^)] 

=-^/^r' A^(z,?)££„(?)exp[-/((0„?+(^„(?))]f/„(z) Jj/exp{-[/((0^^ [VI-33] 
=-Uph-' N{z,t)Y,1^Xt) exp[-/(co„ ?+^„(?))]f/„(z) 2)((0,b-(0„;Yab) 

^ n 

In carrying out this integration we have, again, assumed that changes in £ ^, and (^^^, are 
negligible in a time y~l . Using this result and Equation [ VI- 13a ] we obtain the time 

derivative of the second-order approximation to the time dependence of the population of 
the upper level ~ viz. 



r aa r aa r aa 

= ^.-Y.P..- 



2h 



N{z,i) £££ „ exp {/ [(a)„ -co J?+^„ -(t) J} Uj^z) ^03^^-(0^; y J + c.c. 



1 [ VI-34 ] 



The key point here is that the population has pulsations at the intermode beat 
frequencies! Integrating we find 
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2 

N{z,t)YJl'£.'£r. exp{ i[K-(oJ t+<^,-<^Ji}u^{z)ul{z) 2>(co„-co„; Ya) ^D^b -«m;Yab) + c.c. [ VI-35 ] 

n m 

and since P'bb' =-pl^' with =>Yj, we obtain the second-order approximation to the 
population difference as 

2 

N{z,t)Y.'L'£.'£r. exp{/[(to„-a)J t+^,,-^Ji}u^{z) U*{z) 

n m [VI-36] 

x[®(a)„-to„;Yj + ®(co„-a)„;Y,)][2)(a),,-a)„;Yab) + ®(o>„-a),,;Y,J] 

Substituting this expression into the formal integral of Equation [ VI- 14 ], we find, 
directly, the third-order approximation to the off-diagonal component of the population 
matrix 

2 

N{z,t) IH,^{E^'E.Uiz) K{z) U,{z) 

1 m n 

xexp{-i[(03 -(o„+(oJr+^^-^^+^J}^D((o,b-cOl+03„-(o„;Yab) [VI-37] 

x[2)(a)„-a)„;Ya) + 2t«m-«n;Yb)][^Z^«ab-»n; Yab) + ^^{»m-«ab;Yab)] 

Using Equation [ VI- 10b ] and invoking Equations [ VI-8 ] the general self-consistency 
conditions may be written ^8 



r aa 



in 



K aa K bb 



in 



Pab 



in 



28 It is important to note, that use of Equation [ VI- 10b ] requires evaluation of the following integral : 

L I L 

Numn = \dz N(z) U;{z) mz) Ujz) U„iz) ndz \U,{zt 
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£\ = «k £ k - LLLe i-E „ Im{^kin,„ exp [ / [ VI-38 ] 

1 m n 

%+<^k=^k+<?k-LLL£k'£i£m£„ Re{i^kim„ exp [ / %,^,]] [ VI-39 ] 

1 m n 



with the coefficients summarized in the following table for a standing wave 
configuration: 



Coefficients 


Significance 


«k= 4^ab-»k; Yab) -2 (^JQ^ 


Linear net gain 


<^k= [Kb-«k)/Yab]4«ab-«k;Yab) ^ 


Linear mode pulling 


Pk=^'Kb-»k;Yab) ^3 


Self saturation 


^kl.„ = i^\ r ^1 {l+ [A^.(_) + A^2(»-l) ]/^} ^(«ab -«„; Yab) 


General saturation term 


X [^2^«ab -»„; Yab) + 2^«m -«ab; Yab)] [^X^m -«„; Ya) + Yb)] 





Thus, the coefficients in Equations [ VI-38] and [ VI-39 ] are sensitive to the particulars of the laser mode 
configuration. Since the modal excitations cannot change rapidly , the only integrals of importance are those 
forwhichCO^ -CO; -l-CO^ -a)„ =0 
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t3^ +t3^ = 7 i 

"kkmm^^kmmk * ^ 



Ya+Yb 



pim f {2+ A^,(„_,)/ A^} ^D(a),, -co,; y,,) 



L YaYb 

X{4^ab -«k; Yab) 4^ab Yab) 



+ ■ 



Y, 



ab 



YaYb 



.Ya+Yb, 



[^^«ab k; Yab) + ^^«m -» ab : Yab)] 

X [ 2>(0>m - « k; Ya) + 2^0>m - » k; Yb)] } 



= [(2^'YaYb)/lrf] Im(l^tkmm+1^kmmk) 

= 1^3 {2+ A^2(m-.)/ A^} [4«ab -«k;Yab) ^^ab Yab) 



+ 



2 [Yab] 



YaYb 
Ya+Yb 



Re|^)(a),b k; Yab) [^Wab k; Yab) + ^(Om ab; Yab) 

X [!Z)(a)„ -CO,; Ya) + 2)(co„ -co ,; yj]} 



Stationary coefficient 



Cross-saturation 



%m„ =K-«l+«m-«„) ^+^k-^l+^m-^„ 
Pm = (»ab m)/ Yab ^ ' («ab Yab) ^3 

T^km = [(2^^ Ya Yb)/|pr] +1^ 

kmmk ] 

= 1^3 {2+ A^2(n,-k)/ A^} [(Wab - Wk)/Yab ^ («ab -«k; Yab) ^ (^ab Yab) 



-^[Yabf Im|^D(co,,-co,;Y,,)[25(co,,-co,;Y,,) + ^D(co^-co,,;Y,,)] 

. Ya + Y b 



Relative phase angle 
Self mode pushing 



Cross pushing 
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First-order factor 

and 

Third-order factor 


A^2m= , 1 dzN{z) COS 


\ 1 

z 

1 J 


Spatial factor 


CD^u- v;w) = [i (m- v) + w] 




Lorentzian denominator 


L{u— v; w) = w^ [( ~ + ~] 


-1 


Dimensionless Lorentzian 



Two-Mode Operation; For two modes the amplitude detennining equations ~ i.e. 
Equation [ VI-38 ] - reduce to 

[VI-40] 

2 — 2 1 ^ 2 P 2 "^2 ^ 21 "^ll 



and the frequency determining equations ~ i.e. Equation [ VI-39 ] ~ reduce to 

[ VI-41 ] 



^2+^2 =^2+^2~p2 '^21 A 



W 



Multiplying the amplitude equations by V aY^) obtain the equations of motion 

for the dimensionless intensities 

i, = 2I,[a-^J-%,,I,] 
1^ = 2 I^la^-^^ I2-Q21 A] 
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Stability of possible steady state solutions ; 

Stability criterion for stationary solutions: If /i=//*'+8 and /2=/2^'+e2, ^ and 
82 ^ as t^oo. 

Solution 1: = and 7^'' = a2/p2 so that 

£1 = %][«!- Pi 61-612 (fl2/p2+e2)] 
E2 = 2[a2/p2+e2][«2-p2(«2/p2+e2)-02i ej 



[ VI-43a ] 



or 



82 =-2 (a2/|32) [|32 £2 + 821 eJ+0(8^) 



[ VI-43b ] 



Stability requires that a^= ai-9i2 «2/P2 remain negative. If a{<0 but a^>0 one says 
that I2 inhibits the oscillation of /j (mode inhibition). If becomes large to overcome 
the inhibiting effect of 4 , the solution becomes unstable and /j builds up. 

Solution 2: 9^ and 7^'' 9^ so that 



which has the solution 



a, = 13, 71" + 9,2//' 

«2 = P2 4'' +021 A*'^ 



^W_ (^l/Pl)-(e2l/Pl)(«2/P2) _ ai7Pl 

' 1-C i-c 

^M_ (^2/p2)-(ei2/p2)(^/Pl) _ ^27P2 

' 1-C 1-C 



[ VI-44 ] 



[ VI-45 ] 



R. Victor Jones, March 16, 2000 80 



The Interaction of Radiation and Matter: Semiclassical Theory 



Page 81 



0^, 

where the coupling constant C= . Substituting this solution into Equation 

Pi P2 



[ VI-42 ] we see that 



In matrix forni 



d_ 

dt 



^1 =-2 [pi 81+612 £2] + 0{8') 

£2 = - 2 / [p2 £2 + e 21 81] + 0(8') 



-2 r «; a;(ei2/i3i)i 

l-c[a^e2i/p2) a'^ J 



[ VI-46a ] 



[ VI-46b ] 



The solutions will be stable if the eigenvalues of this equation are both negative. 

^ Det(5^^ -Xl)=0 [ VI-47] 







= X 












£2. 




£2. 



^i + aj 



(■.)■ 



1-C 



1-C 



-4 



^1 ^ 2 

vl^y 



[ VI-48 ] 



Case 1: 



Case 2: 



Weak coupling where C<1 and ai'>0 and a2>0: Stable since 
both eigenvalues negative 

Strong coupling C> 1 and a[<0 and a2<0: Unstable since one 
eigenvalue is positive 
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The van der Pol Negative Resistance Oscillator 



Van der Pol's analysis' of "negative resistance" (e.g., tunnel diode) oscillators prides a valuable 
framework for treating relative simplicity important features of oscillatory systems. 

The characteristic curve of a "negative resistance" device 



LJ 




Voltage 



Consider the following negative resistance oscillatory circuit: 




AVv\— 1 



negative 
resistance 



T 







By simple circuit analysis, it is a straightforword proposition to find the following simple circuit 
equation which is the fundamental van der Pol oscillator equation: 



' B. van der Pol, Radio Rev. 1, 704-754, 1920 and B. van der Pol, Phil. Mag. 3, 65, 1927 
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V (0- [a v(0- P v\t)] V (t)= [ VdP-1 ] 

at a f- 

where tOg = (LC) 

If a is small, it is reasonable to take 

v(t)=^V(t) exp(- i (Op r ) + c.c. [ VdP-2 ] 

Then Equation [ VdP-1 ] becomes without approximation 

j-^ [-tOg V(t)- i 2a)o V{t)+ exp(-/ co,, ? ) + c.c. 

- [a- 3 p v'(r)] 1^ [-/ tOo y(0+ V(0] exp (-/ (A^t) + c.c.j [ VdP-3 ] 

+0)o|i V(Oexp(-/a)o?) + c.c.|= 
If we ignore harmonic generation. Equation [ VdP-3 ] may be approximated as 
- [-12(00^(0+^(0] exp(-jcOor ) + c.c. 



- a 1^ [- j (Oq y( + V{ O] exp (- j cOq ? ) + c.c.j [ VdP-4 ] 

+ ^15 ^|-/(0o| y(0 f H0+ j-jl V(Of y(0]|exp(-/a)oO + t^-t^-j= 
If we make the slow time variation assumption, this equation reduces to 

v(o= ^ [oc- ^ P I y(0 f J v(o= ^ [cc-p I y(o f] v(0 [ vdP-5 ] 

where |3= 3 ]5/4 . This essential Equation [ VI-25b ] in the lecture set entitled The Interaction of 
Radiation and Matter: Semiclassical Theory. We saw there that the general steady state solution 
is given by 

|V(0f = | [VdP-6] 

To study frequency locking we suppose that a driving source (to be precise a current source in 
parallel with the negative resistance) and then the van der Pol equation becomes 
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— J V (0 [a v(0- P v^(0] +03o V (0=0)^ Vq sin o) r 



^— ^ [i exp (-/ 03 r ) + c.c] 



In this case, it is reasonable to take 



1 



v(0= 2^^^-^ exp(-/ CO ? ) + c.c. 
If we again ignore harmonic geneation. Equation [ VdP-7 ] becomes 
[-co' V(0- i 2co V(0+ exp(- i cor) + c.c.j 

- [a- 3 p r)] 1^ [- i co V( 0+ V( O] exp (-/ co r ) + c.c. 
+03 1^ y(0 exp(-j 03r) + c.c.| 

2 

= " [/ exp (- j 03 r ) + c.c] 
Again under the slow time variation assumption, this equation reduces to 

-{o3o-o3')y(0-«coy(0+«co- [a-p I y(of]y(0= ^" 



[ VdP-7 ] 



[ VdP-8 ] 



[ VdP-9 ] 



[ VdP-10 ] 



If we take y(?) = | V[t) | exp(-/ ^[t] ) , this equation separate into the following pair of equations: 



I y(.) 1= i [a-p I y(or] I I- ^ cos 

1 (o3o-03') (0 y 



2 03 



[ VdP-lla] 



[ VdP-llb] 



1 (O3o-03') (0 y 

where d= =(c0o-co) (the "detuning term") and ^=~r777T"i ^^^^ "locking 



2 03 



2|y(.)| 
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.2 a 

coefficient"). For small we can decouple the equations and take =— from Equation 

(0 Vn (0 B 

[ VdP-6 ] so that 1= — | ^ | ~ ~ ^o- If r/^l » 1 relative phase angle changes linearly in 

time at the rate ^ ~ J. As \d/l\ decreases toward unity, the "locking term" subtracts from the 
"detuning term" in one half of a cycle and adds in the other half. At \d/l\ < 1 there are two values 
of the phase angle that yield the "mode locking" condition 4' = — viz. 

I - sin' (d/i) 

'^PL = L. ■ -^( Ji\ [VdP-12] 
[7U+ sm [d/l) 

We can test the stability of these solutions by taking ^{t)= +e {t) and therefore Equation 
[VdP-llb] becomes 

t{t) = /cos (^pl) e{t) [ YdP-13 ] 

and the solutions are stable if 

/cos(^pj<0 [VdP-14a] 

Jl^-d^<0 [VdP-14b] 
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VII. Semiconductor Photonics 

A. PRELIMINARIES; SEMICONDUCTOR BACKGROUND^^ 

THE Crystal Hamiltonian 

For an assembly of atoms the classical energy is the sum of the following: 



• the kinetic energy of the electrons; 

• the potential energy of the electrons in the field of the nuclei; 

• the potential energy of the electrons in one another's electrostatic field; 

• the magnetic energy associated with spin and orbital variables. 

Dividing the electrons into core and valence electrons and leaving out magnetic effects 
leads to the following expression for the crystal Hamiltonian: 



where a and P label the ions, / and m label the electrons, p is the momentum, M is an 
ionic mass, m is the mass of an electron, U (r„ - Rp j is the interionic potential, and 

y (f, - R^j is the valence-electron-ion potential. 



The quantum mechanics of the assembly is treated to a good approximation by taking 
the total wavefunction of the system as the product 



the kinetic energy of the nuclei; 

• the potential energy of the nuclei in one another's electrostatic field; 




[ VII-0] 




[ VII-1 ] 



This discussion draws heavily on B. K. Ridley, Quantum Processes in Semiconductors (3rd edition), Clarendon 
Press (1993). 
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where ^^^r) = 4>(R„,Rp,R5,R^, ...j is the wavefunction of all the ions and 
^(r;R) = ^(r^,,r,,r„,,r„,...;R„,R|3,Rg,R^,...) is the wavefunction of all the 
electrons at the instantaneous ionic positions. The Schrodinger equation is then written 

^(r;R) H 4> (r)+4> (r) H ^(r;R) 

[ VII-2 ] 

+[^ lan^ce ^ (f; R) <J> (r) - ^ (?; R) ^ 4> (r)] = £ ^ (r; r) (r) 
where the total Hamiltonian is parsed into two independent components ~ viz. 

^lauice =L:;^+Lf/(Ra-Rp) . [ vii-3a] 



^.,„ = I^+Iv(f,-Rj+E#|^. [VII-3b] 

I la l,m ri *ml 

The essential assumption of the adiabatic approximation is that the bracketed 
term is negligible and that the global problem may be treated as two independent 
problems ~ viz. 

^lauice ^(r)=^,««,c.^(r) [VII-4a] 

i^,,,,^(f;R)=£,,,,^(f;R) [ VII-4b ] 
As a further refinement, the electron problem must be parsed once more as 

elec ~ {-^elec Astatic {"^ elec } dynamic 
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where ^lec} static defines the problem of the many electron system interacting with the 
static ionic lattice and eiec] dynamic incorporates the effects of the electron-phonon 

interaction. 

LATTICE VIBRATIONS 

For a brief overview of lattice vibrations see the webpage entitled Self Study Materials 
on Solid State Electronics 2 at 

http://deas.harvard.edu/courses/ap216/lectures/ls_2/ls2_u7/sse_tut_2/solid2.html 

Bloch Electrons - Solutions of the Schrodinger equation for a single 
electron or quasi-particle moving in a periodic lattice 

If the electron-electron is averaged out, deviations from the average may be treated 
perturbations — i.e., we make the replacement 

L 1^ _^ I ^ ./..)o + (^ . Ji [ VII-5 ] 

where [pi leads to a constant replusive electronic energy component and {pi 

is a small fluctuating electron-electron interaction. If we neglect these fluctuations, each 
electron interacts independently with a screened lattice potential. In this approximation, 
the electronic wavefunction can be expressed as 

vp(?;R„)= n ¥„{?„; Ro) [VII-6] 

n 

with the requirement that the occupation of one-electron states is in accordance with 
Pauli exclusion principle. Thus, each the wavefunction for each electrons satisfies a 
Schrodinger equation in the form 
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[^^•^)s,a,.c x|/„{f„;Ro)+^(f) x|/„{f„; Ro) = e„ ¥„{?„; Ro)[vii-7a] 

or simplifying the notation 

%^ 

- — V' \|/ (r) + y(r) \|/ (r) = 8 \|/ (r;) [ VII-7b ] 

2 m 



Schematic Periodic Lattice Potential 
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The lattice potential far from the surface of the crystal has the property that, for all 
lattice vectors R = + i 2^2'^ ^ 3 ^3 

y(?)=y(?+R[J [vii-8] 

Bloch's (Floquet's) theorem : The function 

¥k(f)= W"^' Mj^(?)exp(/k ?) [VII-9] 

where "k(^) is a function with the same spatial periodicity as V(f)that is an 
eigenfunction of Equation [ VIT7 ] . 

Proof of theorem: 

First rewrite Equation [ VII-7b ] as 

(r) \|/(r) = 8 \|/(r) [ VII-7b' ] 

Shifting the origin by a lattice vector, we again rewrite Equation [ VITVb ] as 

^(f+R[(]) v(f+R[J=e ¥(f+R[^]) [VII-7b"] 

In light of periodicity expressed in Equation [ VII-8 ], it follows that 

^(?) \|/(f+R[,])=8 \|/(f+R[,]) [VII-10] 

Thus, any linear combination of the possible eigenfunctions 
is a valid eigenstate of the energy £! In particular, let us 
choose the combination 
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= exp( / k r)^\j/(r+Rj,J exp-zk (r+Rj,,J 



[ VII-11 ] 



where k is, for the moment, taken to be an arbitrary complex vector. Since the 

summation includes all possible lattice vectors it must be a periodic function of r 
with the period of the lattice and may be identified with ujr) of Equation 

[ VII- 9 ] .3° If we impose cyclical boundary conditions on '^^{r) , we see that, of 
necessity, k is a real vector! QED 

In light of Bloch's theorem. Equation [ Vll-Vb ] may be rewrite in the form 



^t(f)"-k(f) = 



^(V + /k)Vy(r) 



2m 



"u(f) = £u"-k(f) • [VII-12] 



Since the reciprocal lattice vectors of crystal are defined^i so that exp| ? R[,] ■ j - 1 ^ 
any function written in the form 

exp( /G[„,]-?) 



30 We keep for the record U-J^r) ^ l^^r-l- R[^,]) expf-/ k ■ (r-l- R[^,]) | 

31 The reciprocal lattice vectors: ^[m] =^1^1+ ^2 ^2 '^3 ^3 
where 

^" ^ 3^X3, — 3, X3, — 3,X3- 
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has the spatial periodicity of the lattice. Conversely, we are allowed to expand V{r) 
and Mj,(r) in Fourier series — viz. 

^(f)=L^H exp( /G[,„] ?) [VII-13a] 

[m] 

"t(f)= E"[,„']{k)exp( / G[„,]-?) [ VII-14a] 

[m '] 

where v^^,,] =V' \\\drV (f) exp(-/ G[,„ ,yr) i VII-13b ] 

and %,.i^) = '^' JJJ dfu-^{f)exp(-iG^,„yr) [VII-14b] 

Substituting these expressions into Equation [ VII- 12 ] and equating various Fourier 
components, we obtain the following infinite set of algebraic equations: 

— (G[,„] +k)^ .^[,„]{k)+ i:y[,„.] i^[,„]-[»-]{k) [ VII-15 ] 

which can, in principle, be solved for and M[,„](k) . If we were to turn off the 

lattice potential off to effectuate the empty lattice approximation, the 
eigenfunctions and eigenvalues should, obviously, be those of a free electron and 
given by 

^¥t(f) = M'""' exp ( / k ■ ?) [ VII-16 ] 
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ID EMPTY LATTICE MODEL 



extended zone scheme 



reduced zone scheme 





3tt _2tt 
a a 



a 



TT 

a 



2tt 
a 



3tt 
a 



With gradual restoration of the lattice potential, the wave function is gradually 
transformed from a plane wave, Equation [ VII- 16 ], to a Bloch function. Equation 
[ VII-9 ] . In general, the value of k is well defined since it does not change from its 
original value and the plane wave merely becomes modulated by the function ui.{f). 
Clearly, it is also possible to write the Bloch wave function in the form 



= M~''{ exp ( i G[„] • ?)} exp [ / (k - G[^]) • 



[VII-18] 



and since both u^.{r) and exp( / G[^] r) are spatially periodic, so is their product. 

Thus, the wave function can be considered as obeying the Bloch theorem with any 
empty lattice wave vector k-G^^j. Accordingly, the theorem allows for two 
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alternative modes of classifying solutions ~ i.e. a solution may be specified with the 
original value of k (extended zone scheme) or with that value "reduced" by a reciprocal 
lattice vector plus a band index (reduced zone scheme). 

In the nearly free electron approximation we use the ¥k(r) functions as a basis 
set and take the periodic potential V (r) as a small static perturbation. Using Equations 

[ A-6c ] and [ A-7b ] of this set of lecture notes (at the end of the Section II, Review of 
Basic Quantum Mechanics: Dynamic Behavior of Quantum Systems) we obtain 

1/8- |yl8- \f 
g^^gM ^ e'-e"' [vii-20] 

In light of Equation [ Vll-lSa ] we may write 

\ ^\^l) = i'^)' JJJ exp(-ik'-r)v(f)exp( ik-f) 

= ('^) ' \\\ ^'^P (^'~ ^) ■ ^] ^""P ( ' ■ ^) [ VII-21 ] 

[m] 

= V„5(k + Gj„j-k') 

so that Equations [ Vn-19 ] and [ Vn-20 ] become 

|P_ \~|PP \+V F" \_^ 

l*^k/~|<^k / + ^k+G,„, /pO 



/ £0 Zgo [ VII-22 ] 



k 
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[ VII-23 ] 



[m] ^k+G 



Clearly, the strongest departure from free electron behavior occurs when the 
denominator vanishes — i.e. when 



g g — ^ 

=il 

2m 



+ 2Gr 



[ VII-24 ] 



Equation[ VII-24 ] is the famous Bragg scattering condition or the condition 
which defines the boundaries of the Brillouin zones. At the zone edge for a particular 
pair of degenerate free electron states the appropriate pair of equation from the set in 
Equation [ VII-IQ ] may be approximated 



k+G 



[ VII-25a ] 



or 



8- = - \e° - 



e?^ -e° +4 



k+G 



V, 



[m] 



[ VII-25b ] 



R. Victor Jones, April 4, 2000 



The Interaction of Radiation and Matter: Semiclassical Theory 



Page 92 




R. Victor Jones, April 4, 2000 



The Interaction of Radiation and Matter: Semiclassical Theory 

These considerations then provide the context for the examination of real band 
structures. 

Diamond Structure Brillouin Zone and Nomincature 
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Silicon - band structure 




Density of States Momentum 

10" cm"* eV"' 



For more on band structures go to: 

http://www.deas.harvard.edu/courses/ap216/images/bandgap_engineering/bandgap_engineering.html 

For more on pn junction physics go to: 

http://www.deas.harvard.edu/courses/ap216/images/pnJunction/pnJunction.html 
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B. Optical Properties of Bulk (3-D) Semiconductor^^ 
Free-Carrier Theory 

For most of our considerations a parabolic-band model provides a reasonable 
picture of optical properties. In this model we assume that near a band edge the single 
particle energy of the electron measured with respect to the edge is given by 

^l=-N-=-r^ [VII-26a] 
2m 2m 

where the effective mass m is a measure of the inverse curvature of the band.^^ 
Introducing the notion of a deficiency of electrons in the valence band — i.e. holes -- 
as positive charge carriers, the parabolic-band model is re-interpreted as, respectively, 
the quasiparticle energies of an electron in the conduction band and a hole in the valence 
band 

2m„ 2m„ 



^e,i.=-r-^ = - — [VII-26b] 



8 - = — L-L=^l-!L [VII-26C] 

2m, 2m, 



32 



33 



This section draws heavily on Weng W. Chow, Stephan W. Koch and Murray. Sargent III, Semiconductor- 
Laser Physics, Springer- Verlag (1994) 



That is, the effective mass is given by 



m 



Band edge 



and, thus, in the conduction band the effective mass of electrons is positive and in the 
valence band it is negative — i.e. m^> and m , < 0. 
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where = m^, > and m,^ > accounts for aggregate dynamics of all the electron in 
the valence minus the single empty state. Of course, the single particle picture 
discussed thus far, does not include the critically important effects of the inter-electronic 
Coulomb interactions and, in particular, carrier-carrier scattering . In what is 
usually called free-carrier theory, it is assumed that carrier-carrier scattering 
causes a rapid "thermalization" (relaxation times less than 0. 1 picoseconds) of excited 
conduction band electrons (and valence band holes) and, consequently, it is assumed 
that carriers within a band are in quasiequilibrium with energies distributed 
according to a Fermi-Dirac distribution (see illustration on the next page).^^,, 
Accordingly, the carrier density for a given band {i.e. a= e for the conduction band 
and a= h for the valence band) is determined by the condition 

A^a= ^ E ^ L |expr(3(8„^-^„)j+l| [VII-27] 

states states ^ ^ 

where V is the volume of the sample and |J,„ is the carrier quasichemical 
potential (or imref = quasi-Fermi energy). 3? In the discussion of Equation 



The hole energy is the energy of the filled valence minus the energy of the filled with a vacant electronic state so 
that an increase in hole momentum leads to an increase in hole energy. Since must of necessity include the 
many-body Coulomb interactions among the valence electrons. 

See a Derivation of the Fermi-Dirac distribution function at http://ece-www.colorado.edu/~bart/book/fermi- 
dirac_derivation.htm 

See the applets Fermi Level vs. Carrier Concentration and Doping of Donor and Acceptor Impurities 

(at http://www.acsu.buffalo.edu/~wie/applet/fermi/fermi.html) 
and Fermi Level, Fermi Function and Electron Occupancy of Localized Energy States 

(at http://www.acsu.buffalo.edu/~wie/applet/fermi/functionAndStates/functionAndStates.html). 

As are the electron and hole energies ^ , the imrefs are measured with respect to the appropriate band 
edge — viz. 

\^.={Ef)c-Ec and [l, = Ey-{E^)^ 
Thus, if < the band does contain enough carriers to populate any state with a probability greater than one-half. 
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[ IV-5 ] of this set of lectures, it was argued that the 3-D density of states in k space is 
given by 



, , V 2 

p[k) dk= — k dk 



[ VII-28 ] 



which translates (within a given band) into a density of states per unit energy per unit 
volume of 



dk 1 



Y ^ ' dE . 2r 

ak 



ak 



[ VII-29] 



Hence, for a 3D semiconductor. Equation [ VII-27 ] becomes 



N^= Jp(8„)/(e„)je„=^fef Jy8:{exp[p(8„-^J}"'j8„ [VII-30] 



Fermi FuncliDn 




Er-nj.nv 



£f- D.1iV 



1(E) 1 JO 0.5 




- Ef + D.ltV 



Ef - Cl.l 



1. 




HE) I.D 0.5 
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Optical Matrix Elements 

Since optical interactions in a semiconductor are essential distributed, it is probably 
more precise and appropriate in this case to write the interaction Hamiltonian as ^8 

i?/- J =- A • !P = - A • => ^ [e • !P ] exp f i (co,,, t- k,,,- ?)] + . . . [ VII-31 ] 

m m 2m >■ ^ 

Thus, the critical optical matrix element ^9 between a state in the valence band and one 
in the conduction band is 

{^4 ^ I / = exp[ i(0,^,t] M^- ^ [ VII-32a ] 

where the so called frans/f/on matrix element is given by 



V He k 



,^\df /-(?) exp[-^ k ■ ?] {e-T" exp[-/ k„,- ?]} u -,(?) exp[ / k'-?] [ VII-32b ] 
Using the expansion for in Equation [ Vn-12a ] the matrix element becomes 
Mvk, ok' = L L"v[m](k) Wc[«-](k') 

[ml [ml 

^ ' [VII-33] 
X \\\df exp[-/ (k+ k,,, + G[„,)- r] {e • ^} exp[ i (k'+ G^^,, ) • r] 



38 In the next set of lectures, The Interaction of Radiation and Matter: Quantum Theory, we develop the 
following expression for the nonrelativistic Hamiltonian of a single charged particle: 

^ = ^"^^^ -?A ]' 
where fP is the canonical conjugate momentum of the charged particle. 

39 The Fermi Golden Rule and first-order perturbation theory teUs us that the radiation-induced transition 
rate is given by 

If l</K.JO|s(e/-e,)- 
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However, !P ^-ihV so that 



[m] [m'] 



X 



III J r exp [-/ (k + k,,,- k'+ G[^] - G[„,]) • r 



[ VII-34 ] 



With the parsing r = R^fj + p and using the defining relationship exp( / Rj^j ■ G[,„]) =1 , 
the volume integral may be expressed 



d r {r} exp [ / (k- k - k^^,+ G[,„ - G[„,]) ■ r 

= lll^f {P}exp /(k'-k-k,,,+ G[,„,]-G[,„]) p 

xj^exp[ /{k'-k-kj R[^ 



unit cell 



[ VII-35 ] 



The summation I^exp^ i (k'- k - k^.^^) ■ R[ 



unless 



k'= k + k^^,-k. 



[ VII-36 ] 



i.e. essentially only vertical transitions are allowed! Therefore, 



M - - ^5-- 

V He k kk 



y. 



sample 



V^unit cell J [m] [m ' 



Il{^^(k + G,„J}<f^,{k)^.^[„,.j{k) 
X III Jr exp [/(G[„, .J -G[,J.p 

unit cell 



.[ VII-37a] 



For a given polarization of the applied field, the transition matrix can be written 
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He i' = {polarization factor} M . [ VII-37b ] 

where the polarization factors are of order one and typical values of |M| are given in 
the following table:'*'^ 



Material system 


2|M|7m (ineV) 


Ga As 


28.8 


AlxGai.xAs (x<0.3) 


29.83 + 2.85 x 


Inx Gai-x As 


28.8 - 6.6 X 


InP 


19.7 


Ini.x Gax Asy Pi.y (x = 0.47 y) 


19.7 + 5.6 y 



Of course, the radiation must also satisfy the energy conservation condition 
AEc = h(Oi =8„ + + =£„ +■ 



where 



2m ^ 2n\ 
J_ _ J_ J_ 

m, m„ niu 



2m, 



[ VII-38a ] 
[ VII-38b ] 



"FREE-CARRIER THEORY" OF OPTICAL PROCESSESS 

The bottom line is that in the free-carrier theory of optical interactions, the 
effective Hamiltonian for the carriers ~ i.e. "free" particle kinetic energy plus 
electromagnetic interaction ~ is separable into a series of k -dependent terms.^^i Thus, 



From Diode Lasers and Photonic Integrated Circuits by Larry A. Coldren and Scott W. Corzine, Wiley (1995) 
In particular, it can be shown - see Weng W. Chow, Stephan W. Koch and Murray. Sargent HI, 
Semiconductor-Laser Physics, Springer- Verlag (1994) — that 
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if we neglect kk' correlations in the treatment of optical properties, we need only 
consider the following manifold of states: 





1,1^: t) 









2)E|l^0^:i) 




|0)E|o,0,:i) 



l)E|o,l,:h) 



Equation [ VII-37a ] tells us that the absorption or emission of a single photon connect 
only the two states in the manifold which can be identified as |0,, 0^, : k^ and |l,, l,^ : k/ 

where the first number in the ket specifies electron occupancy of a state with 
momentum ti k and the second hole occupancy of a state with momentum - ^ k . 



2m* 



2m 



hv 



^-k ^-k - [^^k 4 ^k + «k ^-k ] E(z, r) ( 



where {fl^, flj} and{^j,^g} are, respectively, electron and hole {creation, destruction} operators. |Ij is 
the dipole matrix element between vertical states in the valence and conduction bands. 
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ENERGY DIAGRAM FOR BULK LASERS 



Electron Energy 




{exp[p(Eh- Ef^)] + 1} 



Since no kk' correlations are involved in the free carrier theory, the density operator of 
the complete system may be expressed as a product of component density operators 

P=np(k) [VII-39] 

u 

so that the Schrodinger equation of motion for the density operator is also separable — viz. 
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p(k) = 4[^-(k)'pW]- 



[ VII-39 ] 



C. Injection Laser Theory 

(For pictorials on laser configurations go to: 

http://www.deas.harvard.edu/courses/ap216/images/injection_lasers/laser_structures/laser_structures.html) 

A fairly satisfactory model of lasing (and other optical processes) in semiconductors 
may be obtain by adaptation of the two-level, semiclassical discussed earlier. In light 
of the discussion in the previous section, we adapt the two-level theory by making the 
identifications 

pjk) , ,{k) = ( 1,. 1, :k I p I 1, l,:k ) = ( 3:k I p I 3:k ) = P33 (k) 

P4k) ^Poooo{k) = (0,0,:k |p|0, 0,:k)=(0:k |p| 0:k ) = Po„ (k) [VII-40] 

p4»^) ^Pnoo{k) = (l.l„:k |p|0, 0,:k)=(3:k |p|0:k)==P3o (k) 

Given the free carrier effective Hamiltonian discussed above and the equation of 
motion expressed in Equation [ VII-31 ], we can write 

P3o(k) = (3:k|p{k)|0:k) 

= -^[( 3:k I ^,,,{k) |3:k)-( 0:k | i/f,,,{k) | 0:k )] p,, (k) 

+ i(3:k |i/f4k)|0:k)[ p33{k)-p„o{k)] 

P33 (k) = ( 3:k I p{k) I 3:k ) = i [ P3,, (k) ( 0:k | ^ ,,,{k) | 3:k )] 

-^[P03{k)(3:k|i^e«{k)|0:k)] 



[ VII-41a ] 



[ VII-41b ] 
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Poo {k)= ( 0:k |p(k) I 0:k ) = - i [ (k) ( 0:k I ^ ^^k) | 3:k )] 

^ [ vn-4ic ] 

+ ^[po3{k)(3:k|i^,^(k)|0:k)] 

which tells us, once again, that p^^ ^^pj^ is driven by p^^ -p^^, =>p33 -Poo ^^^^ 
versa. There is, however, an important new element to be considered in this problem to 
wit, as we have illustrated above, we are, in fact, actually dealing with a four-level 
system which includes of the states 1 1^ 0;,:k ^ and | 0^ l^:k ^ as well as the pseudo "a" 

and "b" states. Fortunately, inclusion of these states does not unduly complicate the 

analysis since, as discussed above, electron-electron scattering induces a rapid relaxation or 

"thermalization" of the probability of finding a given state occupied. After thermaUzation, 

the probability of finding a state with a particular momentum value in a given band is 

provided by a Fermi-Dirac distribution referenced to a quasi-Fermi energy or imref 
appropriate to that band. To obtain an expression for p aa ~P bb we note ihaAP- 

tr[p(k)]=Poo(k)+Pij(k)+p,2(k)+p33(k)=l [ vn-42] 

so that 

Paa-Pbb=p33 (k) - P oo (k) = { P33 (k) +P22 (k) } + { P33 (k) +p,i(k) } " 1 
= n,(k) + n, {k)-l 



[ Vll-43 ] 



where n^(k) is the probability'^^ of finding a electron with momentum h k independent of 
whether or not there is a hole with momentum - k and n^, (k) is the corresponding 



That is, there must be either or 1 electrons and or 1 holes in the given k state. 

In the second quantization representation of the problem, the n' s are expectation values of the appropriate 
number operators. 
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probability for a hole. Adapting Equations [ VI-13a-c ] from earlier in this set of lecture 
notes, we may write 

< (k) = A,-,- Y„r (k) - T- (k) n, (k) - {k)[ ^ - [ i 1/^, (z,t; k) p (k) + c.c] [ VII-44a ] 

< (k) = A,-,- Y„r n, (k) - r- (k) n, (k) - n, (k)[ ^ - [/ 1/^, (z,t; k) p (k) + c.c] [ VII-44b ] 
P ab{k) =-(i CO,,, +Yab) P ab (^) + * '^abl^'^ k) (k) + (k) " l] [ VII-44c ] 

where A^- [ a= e or h) represents the pumping rate due to carrier injection and n„ (k)| 

the effect of carrier-carrier scattering. Y„n^ (k) and I^n^(k)n^(kj are, respectively, 

phenomenological representations of nonradiative decay and radiative "bimoleculai^' 
recombination (spontaneous emission), respectively. The probability difference (gain 
factor) 



A(k)^n^(k) + n,(k)-l = n^(k)n, {k)|l- 



f r 1 T 1 11 



-1 



-1 



}■ [ vn-45 ] 

J 



is the critical factor in the analysis of stimulated processes - i.e. for inversion 
A(k)>0. This gain factor may be written as A(k) = /i^(k)-[l -n^(k)j ~ i.e. the 

population inversion is proportional to the probability difference of an electron in 
corresponding k -states in the conduction and valence bands. Alternatively and more 
usefully, we see that it varies directly with the product of 



44 That is 



Paa = ^a- Ya Paa " ^"'^^ab P.^+^c] 
Pbb = ^b - y„ Tbb + [i '^ab Tab + C.c] 
Pab = -( i + Yab) P ab + ^' '^^ab [paa ' f^bb ] 



R. Victor Jones, April 4, 2000 



The Interaction of Radiation and Matter: Semiclassical Theory 



PAGE 106 



the spontaneous emission factor 

and the so called absorption factor 



1- 



«.{k)«/,(k) 

r 1 .T 1 1. 



-1 



-1 



In the spirit of the free carrier model and in the unsaturated limit the ^( (k) 's are given 
by the quasi equilibrium Fermi-Dirac distributions — viz. 



"4k)^/„.= fxp[|3(e„,-^l„) j+1 



[ VII-46 ] 



Thus, the absorption factor is given by 

r 1 .T , 1 



1- 



-1 



(k) n, (k) 



-1 



l-exi[|3(8^,+8,,-|a,-^,)j 
l-exp[|3{8,-A^)] 



[ vn-47 ] 



and the spontaneous emission factor is given by 

(k) n, (k) = I expf-P jexp [ P{8-,- A^l)] +1 } |exi[|3(8,;^- ^,)j+l} [ vn-48 ] 

where8,= ^ and A|I= ^1, + ^,= (£^)^-(£^)^-8, . « 



For plotting purposes it is useful to write the complete gain factor as 







jexp 


V 



X r 

h " Xp ~\~ X^ 



+ 1 \ -iexp 



vl + r 1 + r 

where r= mjm^, x= fKO-jE^, x,= A£^^ ,and x, = \i,jE 



X-\ hX 



1+r 1+r 



+ 1 
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absorption spontaneous 




-2 -1 1 2 3 4 5 



gain 

P(P'h- ^h) values 2.0 




Therefore, for population inversion A|J,— 8 — {Ef)^- {Ep)^ - (£^+£) > which is 
a very stringent condition! 
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Adapting Equation [ VI-26a ] from earlier in this set of lecture notes,^^ we may see that 
the small signal gain in a semiconductor, may be expressed in the form 

Gain {(£>-^ \^y'k,c k| ^(^) density of paired statesat k} X^O)- -0) y^,,) [ Vn-49 ] 

in buLc material the density of paired states varies as the square root of the energy. 

Therefore, we may draw the following gain curves: 



. 14 




Gain in a bulk semiconductor with A / E = 1 .00, 1 .02, 1 .04, 1 .06, 1 .08, and 1 . 10 



46 That is 



1 
2 
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GAIN VARIATION WITH CARRIER CONCENTRATION 



1.5 



c 0.5 



01 

> 

0) 



-0.5 



increasing 
injection 
current 




Quantum Confinement 

In an idealized quantum well V[z) = 

which yields a single electron energy 



[0 \z\<Lj2 
[oo \z\>Lj2 



^ /- \ 7c' 



2m^ 4 



2 m, 



and density of paired states which is a constant or a staircase. 
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0.98 1.02 1.04 1.06 1.08 1.1 1.12 



Gain in a 2D semiconductor with A Ep / E = 1 .06, 1 .08, and 1.10 



O.S r 
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